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INTRODUCTION: 


Magnetic  Resonance  Imaging  (MRI)  is  an  excellent  imaging  modality  for  many 
conditions,  but  to  date  there  has  been  limited  success  in  harnessing  this  modality  for  the  guidance 
of  interventional  procedures.  MRI  is  an  ideal  interventional  guidance  modality:  it  provides  near 
real-time  high-resolution  images  at  arbitrary  orientations  and  is  able  to  monitor  therapeutic 
agents,  surgical  tools,  biomechanical  tissue  properties,  and  physiological  function.  At  the  same 
time,  MRI  poses  formidable  engineering  challenges  by  severely  limited  access  to  the  patient  and 
high  magnetic  field  that  prevents  the  use  of  conventional  materials  and  electronic  equipment. 
Currently,  no  technological  solution  exists  to  assist  MRI  guided  prostate  interventions  in  an 
accurate,  simple,  and  economical  manner.  The  objective  of  our  research  is  to  make  conventional 
diagnostic  closed  high-field  MRI  scanners  available  for  guiding  prostatic  needle  placement 
interventions.  Our  approach  is  to  employ  robotic  assistant  devices  for  assisting  clinicians  with 
needle  placement  while  taking  advantage  of  real  time  MR  images  for  guidance.  In  this  work,  we 
apply  a  remotely  actuated  robotic  needle  placement  system  in  standard  closed  high-field  magnets, 
optimized  for  transperineal  prostate  biopsy  and  brachytherapy  seed  placement. 
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BODY: 


Please  see  the  publications  included  in  the  appendix  for  a  detailed  account  of  the 
published  work  to  date.  Specific  accomplishments  in  Year  2  by  aim  of  the  research  plan  are 
detailed  below. 

Aim  1:  Refined  Requirements  and  System  Design 

In  Year  1,  we  worked  closely  with  our  collaborators  at  the  Brigham  and  Women’s 
Hospital  and  Johns  Hopkins  University  to  refine  the  requirements  of  the  robotic  system  [Song 
2010],  The 

most  critical  parameters  relevant  to  design  of  the  robotic  assistant  are  related  to  patient 
positioning  and  placement  of  the  robot. 

In  Year  2  we  focused  on  the  workflow.  To  satisfy  the  design  requirements  of 
interventional  needle  placement  robot  we  proposed  a  system  workflow  that  mimics  traditional 
transrectal  ultrasound  (TRUS)  guided  prostate  needle  insertions  to  maximally  retain  surgical 
workflow  that  physicians  are  comfortable  with.  This  workflow  as  shown  in  Fig.  1  below  includes 
five  states  of  operation  and  follows  a  coherent  procedure. 

•  Initialization.  The  hardware  and  software  system  is  initialized.  In  this  state,  the  operator 
prepares  the  robot  by  connecting  the  robot  controller,  applying  sterile  drape,  attaching 
sterile  needle  contacting  components,  and  inserting  the  needle.  The  robot  is  calibrated  to  a 
pre-defined  home  position  and  loads  the  robot  configuration  from  XML  file. 

•  Planning.  Pre-operative  MR  images  are  loaded  into  the  3D  Sheer  and  the  target  locations 
are  selected  in  patient/image  RAS  coordinates. 

•  Calibration.  A  series  of  transverse  images  of  a  fiducial  frame  on  the  robot  are  acquired. 
Multiple  images  are  used  to  perform  multi-slice  registration  to  enhance  system  accuracy. 
This  registration  step  enables  conversion  between  image  coordinates  to  robot  coordinates. 

•  Targeting.  Needle  target  is  selected  from  the  Sheer  software  and  this  desired  position  is 
transmitted  to  the  robot  interface  software  module  to  process  inverse  kinematics  and  the 
calculated  joint  command  is  used  to  drive  piezoelectric  motors.  Targets  may  also  be 
directly  entered  or  adjusted  in  the  robot  proxy  interface.  Real-time  MR  images  can  be 
acquired  during  insertion  that  enable  visualization  of  the  tool  path. 

•  Verification.  The  robot  forward  kinematics  calculates  actual  needle  tip  position  (from 
encoder  measurements  and  registration  results)  which  is  displayed  in  the  3D  Sheer.  Post 
insertion  MR  images  are  acquired  and  displayed  with  overlaid  target  and  actual  robot 
position. 
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Figure  1  -  Workflow  of  MRl-guided  needle  placement  that  mimics  that  of  traditional  TRUS- 

guided  prostate  needle  intervention 

Aim  2:  Advanced  Needle  Driver  Development 

Actuation 

In  Year  1  we  developed  a  prototype  piezoelectric  actuation  system.  By  leveraging  our 
previous  effort  on  piezoelectric  actuation  and  motor  driver  development,  in  Year  2  we  further 
refined  the  piezoelectric  driver  design  as  shown  in  Fig.  2  [Cole  2001,  Cole  2011,  Fischer  2011], 
There  are  many  shapes  and  forms  of  piezoelectric  actuators,  and  they  generally  fall  into  two 
major  categories:  harmonic  and  non-harmonic.  While  non-harmonic  motors,  operate  at  a  much 
lower  frequency  than  harmonic  motors  (750//z  to  3k Hz)  they  require  a  high  precision  non- 
sinusoidal  waveform  to  be  operated  most  effectively.  Although  we  have  focused  on  PiezoMotor 
nonharmonic  actuators,  we  have  shown  our  drive  system  to  operate  with  all  popular 
commercially  available  piezo  actuators.  The  signal  processing  stage  receives  an  input  in  the  form 
of  a  velocity  or  position  set  point.  This  signal  processing  block  can  then  create  four  independent 
channels  of  driving  waveforms  with  a  data  rate  of  12.5  mega  samples  per  second  of  14  bit 
resolution.  The  signal  generators  include  a  high  power  output  amplification  stage,  which  passes 
its  signals  out  to  the  actuators  through  a  7i  filter.  The  boards  also  integrate  position  encoding 
(from  differential  encoder  receivers)  such  that  closed-loop  motion  control  loops  can  take  place 
within  the  FPGA,  distributing  the  real-time  processing  of  the  system.  Fig.  2  shows  the  latest 
board  design  and  the  backplane  configuration  that  goes  inside  the  controller  enclosure  containing 
five  piezoelectric  driver  boards. 
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Figure  2  -  Latest  iteration  (Rev.  C)  of  the  piezoelectric  actuator  driver  architecture.  This  custom 
drive  circuit  provides  four  FPGA-based  waveform  synthesizers  and  corresponding  amplifiers 
that  read  motor-specific  waveforms  from  an  SD  card(left).  Multiple  driver  boards  can  be  utilized 
on  a  single  backplane  for  easy  connection  to  the  robot  and  enables  communication  through  a 

single  fiber  optic  Ethernet  connection  (right). 


Sensing 

Sensing  is  required  to  implement  a  teleoperation  framework  with  force  feedback  [Su 
2011c].  In  Year  1  we  developed  a  3  DOF  fiber  optic  force/torque  sensor  for  measuring  tissue- 
needle  interactions  [Su  2009,  Su  2010].  In  Year  2,  we  further  designed  a  Fabry-Perot 
interferometry  (FPI)  based  miniaturized  force  sensor.  Even  though  fiber  optic  sensors  are 
generally  accredited  for  better  survivability  in  hazardous  environments,  meticulous  design 
considerations  are  still  required  to  build  robust  and  durable  sensors.  For  intra-body  applications, 
sensor  miniaturization  and  tool  integrability  are  essential  for  surgical  procedures.  In  particular, 
for  prostate  interventions,  the  sensor  should  be  small  in  size  and  easy  to  integrate  with  surgical 
needle. 


Figure  3  -  Fabry-Perot  sensing  principle:  light  propagation  in  Fabry-Perot  cavity  (left)  and 

resulting  fringe  pattern  (right). 


Figure  4  -  Magnified  fiber  FPI  strain  sensor  with  three  segment  dimensions  (left)  and  FPI 
embedded  in  an  ABS  cantilever  beam  and  the  inset  shows  the  fiber  with  a  cent  (right). 
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Besides  immune  to  electromagnetic  and  RF  signal  and  substantially  cheaper  than  fiber 
Bragg  grating  sensor,  the  advantages  of  FPI  sensor  includes:  1)  static/dynamic  response 
capability,  2)  high  sensitivity  and  resolution,  3)  no  interference  due  to  cable  bending  and  4) 
robust  to  a  large  range  of  temperature  variation  (-40°~250°)  due  to  air  gap  insulation  to  the 
sensing  region.  As  shown  in  Fig.  3  &  4,  the  main  component  of  the  FPI  is  the  sensing  cavity, 
measuring  15.8pm  wide.  A  glass  capillary  covering  the  sensing  region  is  fusion  welded  to  the 
fiber  in  two  locations  and  encapsulates  the  sensor.  There  is  an  air  gap  of  approximately  100.5  pm 
wide.  The  total  length  of  the  FPI  sensor,  including  the  glass  capillary,  and  bare  fiber  is 
approximately  20  mm. 

Fig.  5  (top)  shows  the  configuration  of  the  custom-made  FPI  sensor  interface  and  sensor 
location  on  the  robot  needle  driver  [Su  201  Id].  Due  to  negligible  friction  force  between  the 
needle  and  needle  guide,  the  reaction  forces  between  the  mechanism  (top  plate)  and  the  actuator 
drive  rod  is  used  to  measure  needle  insertion  force  as  shown  in  the  Fig.  5  (bottom-right).  The 
beam  to  hold  the  actuator  rod  has  a  small  1:59mm  groove  to  lay  the  sensor  that  would  extend 
30mm  along  its  side  where  the  FPI  sensor  could  be  embedded.  10  Newton  interaction  force  is  the 
maximum  required  for  each  sensor.  The  finite  element  analysis  illustrates  the  maximum  strain 
lOOps  under  10  Newton  axial  force  using  ABS  plastic  material.  Calibration  results  are  shown  in 
Fig.  5  (bottom-left). 


Figure  5  -  Fabry-Perot  interferometry-based  force  sensing  of  needle  insertion  force.  Intended 
configuration  (top),  calibration  of  a  test  fixture  (bottom-left)  and  FEM  of  the  needle  driver  to 
predict  deflection  under  typical  10N  insertion  force. 
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Needle  Manipulation 

We  have  developed  the  first  prototype  of  a  MRI-compatible  piezoelectric  actuated  robot 
integrated  with  a  high-resolution  fiber  optic  sensor  for  prostate  brachytherapy  with  real-time  in 
situ  needle  steering  capability  in  3T  MRI  [Su  2011c].  The  6  DOF  robot  consists  of  a  modular  3 
DOF  needle  driver  with  fiducial  tracking  frame  and  a  3  DOF  actuated  Cartesian  stage.  The 
needle  driver  provides  needle  cannula  rotation  and  translation  (2  DOF)  and  stylet  translation  (1 
DOF).  The  driver  mimics  the  manual  physician  gesture  by  two  point  grasping  (stylet  and 
cannula).  The  robot  design  is  shown  in  Fig.  6. 


Fiducial  x7  Optical  encoder 


Fidu< 


Figure  6  -  Physical  prototype  of  6-DOF piezoelectric  actuated  needle  placement  robot  (left).  An 
exploded  view  of  the  needle  clamping  mechanism,  fiducial  frame  and  motor  fixture  (right). 

Aim  3:  Robotic  System  Integration 

We  have  developed  an  integrated  system  consisting  of  four  major  components: 
Navigation  software,  robot  interface  software  and  communication  protocols,  robot  controller,  and 
the  robot  itself  [Su  2011a],  Three-dimensional  surgical  navigation  software  3D  Sheer  serves  as  a 
user  interface  to  visualize  planning  MRI  images  and  define  targets  in  image/patient  RAS 
coordinates.  We  have  developed  a  robot  interface  software  that  replies  upon  OpenIGTLink,  a 
network  protocol  designed  to  work  on  the  TCP/IP  application  layer  to  handle  transform,  image 
and  device  control.  Communication  from  the  robot  interface  to  the  robot  controller  is  through  a 
fiber  optic  Ethernet  connection  run  through  the  patch  panel  wave  guide,  as  this  eliminates  a  large 
source  of  noise  that  is  introduced  when  electrical  signals  are  passed  through  the  walls  of  the 
scanner  room. 

Fig.  7  shows  the  system  architecture  -  Sheer  is  used  to  select  the  target  position  in  MR 
images,  the  target  is  transmitted  through  OpenIGTLink  to  the  robot  interface  software, 
registration  results  and  inverse  kinematics  are  used  to  determine  the  appropriate  robot 
configuration,  the  commands  are  sent  to  the  robot  controller,  and  the  controller  actuates  the 
piezoelectric  motors  under  precise  closed  loop  control.  For  verification,  the  tip  position  is 
calculated  based  on  the  robot  forward  kinematics  and  registration  results  to  overlay  the  actual 
needle  location  (as  determined  from  robot  encoders)  on  the  MR  images  in  Sheer. 
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The  complete  robotic  system  has  been  tested  in  a  3T  MRI  scanner  during  in  phantom 
trials.  The  robot  shown  in  Fig.  8  (left)  resides  on  the  scanner  bed  and  operates  within  the  bore  of 
the  scanner  without  any  noticeable  MR  image  degradation.  The  robot  is  connected  via  a  single 
shielded  electrical  cable  to  the  robot  controller  which  sits  in  the  scanner  room  beside  the  bed  as 
shown  in  Fig.  8  (right). 


Navigation  Software 
(3D  Slicer) 


I  target 

I _ EffisMgn 


Control  PC 

□z 


QpenlGTUnfc 


r  Registration 


Piezoelectric 

Motor 


Figure  7  -  System  architecture:  system  architecture  diagram  shows  the  connection  and  data 
flow.  OpenIGTLink  is  used  to  exchange  control,  position,  and  image  data  (left).  The 

corresponding  components  (right). 


Figure  8  -  A  detailed  view  of  the  robot  with  phantom  inside  MRI  (left).  The  prototype  needle 
placement  robot  for  percutaneous  prostate  interventions  inside  3T  MRI  scanner  and  the  robot 
controller  is  inside  the  scanner  room  beside  the  scanner  (right). 
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Aim  4:  Validation 


MRI  compatibility  of  the  initial  prototype  was  evaluated  in  Year  1.  We  have  continued  to 
confirm  compatibility  of  the  latest  controller  design,  and  found  that  even  during  simultaneous 
robot  motion  and  imaging,  no  more  than  5%  SNR  loss  is  found.  Therefore,  we  are  confident  in 
the  ability  of  the  current  system  to  operate  during  imaging  with  no  noticeable  image  quality 
degradation. 

With  the  fully  functional  system,  a  series  of  experiments  are  performed  to  evaluate  the 
system  performance  for  needle  insertion  and  steering  capability  under  real-time  3T  MRI- 
guidance  [Su  2011c].  The  first  test  is  the  dynamic  needle  insertion.  Ballistics  gelatin  is  utilized  as 
a  tissue  phantom  for  evaluating  needle  insertion.  This  imaging  protocol  provides  2.5Hz  update 
rates.  Needle  insertion  precisely  controlled  by  the  robot  controller  during  simultaneous  imaging. 
Fig.  9  (left)  depicts  six  needle  insertion  snapshots  during  3T  echo-planar  imaging  at  0.4  second 
interval.  The  needle  shaft  and  tip  trajectories  are  clearly  visualized  in  the  phantom  image  without 
major  interference  during  robot  motion. 


Figure  9  -  Needle  insertion  snapshots  during  real-time  3T  echo-planar  imaging  at  0.4  second 
interval  (left)  and  Bevel  tip  needle  steering  image  in  3T  MRI  (right). 


In  the  second  test,  the  same  22  Gauge  medical  needle  is  used  to  demonstrate 
compensation  of  the  needle  path  by  controlled  bevel-tip  steering  with  MRI  visualization.  The 
bevel  tip  is  rotated  toward  left  before  insertion.  T2-weighted  fast  spin  echo  (field  of  view 
240mm,  echo  time  90ms,  repetition  time  3000ms,  flip  angle  90°)  illustrates  the  final  needle  shape 
and  tip  position.  The  same  procedure  is  repeated  for  bevel  right  before  insertion  and  the  results 
are  shown  in  Fig.  9  (right).  There  is  no  visually  identifiable  interference  during  needle  robot 
controlled  insertion. 


To  determine  the  systems  effectiveness  at  controlling  actuators,  it  was  designed  to  drive  a 
PiezoLEGS  rotary  motor  that  was  coupled  to  a  US  Digital,  1250  count  quadrature  optical 
encoder,  and  then  driven  through  a  series  of  step  response  and  sine  response  experiments.  These 
experiments  were  structured  such  that  a  position  set  point  was  supplied  to  the  driver  on  a  timed 
schedule,  and  the  encoder  information  was  read  and  recorded  in  real  time  as  shown  in  the  figure 
below.  These  results  shown  in  Fig.  9  demonstrate  the  ability  of  the  controller  system  to  drive  an 
actuator  with  finer  precision  than  a  single  encoder  tick,  which  in  this  case  corresponds  to  a  0.072 
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degree  or  change.  Similar  results  or  encoder-level  accuracy  were  found  with  the  linear  motors 
which  are  measured  using  a  1000  cpi  quadrature  encoder,  thus  an  accuracy  of  0.006mm.  Note 
that  the  purpose  of  this  study  was  to  measure  the  inherent  accuracy  of  the  robot  actuators,  the 
limiting  factor  in  the  clinical  system  is  the  pixel  size  of  the  selected  MRI  scan  protocol  and  the 
corresponding  registration  accuracy  (typically  sub-pixel).  Phantom  studies  are  run  with  0.5mm 
pixel  size. 


Figure  10  -  Motor  position  (blue)  and  set  point  (red)  plotted  against  time  for  two  trials:  dynamic 

sinusoid  motion  (left)  and  step  response  (right). 
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KEY  RESEARCH  ACCOMPLISHMENTS: 


•  Defined  system  workflow 

•  Developed  next  generation,  clinically  viable  control  system 

•  Developed  an  integrated  networked  hardware  and  software  system  to  support  the 
proposed  surgical  workflow 

•  Refined  6  DOF  needle  placement  robot  including  3  DOF  needle  driver  module 

•  Developed  a  miniaturized  Fabry-Perot  interferometry  based  fiberoptic  force  sensor 

•  Demonstrated  a  fully  functionally  robot  inside  3T  MRI  scanner 

•  Validated  capability  of  needle  placement  under  real-time  MR  imaging 

•  Evaluated  robotic  system  accuracy 

REPORTABLE  OUTCOMES: 

Please  see  References  section  for  publications  resulting  from  this  research. 


CONCLUSION: 

Based  on  the  requirements  defined  in  Year  1,  we  defined  a  surgical  workflow  and 
developed  an  integrated  hardware  and  software  system  that  consists  of  a  refined  6  DOF  needle 
placement  robot  with  fiber  optic  force  sensing  utilizing  a  miniaturized  Fabry-Perot  interferometry 
sensor,  a  software  controller  with  Ethernet  communication,  and  an  improved  piezoelectric  motor 
controller.  With  ex  vivo  phantom  study,  we  have  integrated  the  system  and  evaluated  system 
workflow,  MRI  compatibility,  image  registration  accuracy  and  robot  tracking  accuracy.  In  the 
coming  year,  we  will  enhance  system  accuracy  and  reliability,  complete  bilateral  teleoperation 
with  fiberoptic  sensor  based  force  feedback.  The  final  deliverable  at  the  end  of  the  project  period 
will  be  a  teleoperated  version  of  the  robot  shown  that  is  tested  in  phantom  studies  in  preparation 
for  pre-clinical  trials  at  BWH. 
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for  the  proposed  robotic  system. 
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Abstract —  This  paper  presents  the  first  prototype  of  a 
magnetic  resonance  imaging  (MRI)  compatible  piezoelectric 
actuated  robot  integrated  with  a  high-resolution  fiber  optic 
sensor  for  prostate  brachytherapy  with  real-time  in  situ  needle 
steering  capability  in  3T  MRI.  The  6-degrees-of-freedom  (DOF) 
robot  consists  of  a  modular  3-DOF  needle  driver  with  fiducial 
tracking  frame  and  a  3-DOF  actuated  Cartesian  stage.  The 
needle  driver  provides  needle  cannula  rotation  and  translation 
(2-DOF)  and  stylet  translation  (1-DOF).  The  driver  mimics  the 
manual  physician  gesture  by  two  point  grasping.  To  render 
proprioception  associated  with  prostate  interventions,  a  Fabry- 
Perot  interferometer  based  fiber  optic  strain  sensor  is  designed 
to  provide  high-resolution  axial  needle  insertion  force  measure¬ 
ment  and  is  robust  to  large  range  of  temperature  variation.  The 
paper  explains  the  robot  mechanism,  controller  design,  optical 
modeling  and  opto-mechanical  design  of  the  force  sensor.  MRI 
compatibility  of  the  robot  is  evaluated  under  3T  MRI  using 
standard  prostate  imaging  sequences  and  average  signal  noise 
ratio  (SNR)  loss  is  limited  to  2%  during  actuator  motion.  A 
dynamic  needle  insertion  is  performed  and  bevel  tip  needle 
steering  capability  is  demonstrated  under  continuous  real-time 
MRI  guidance,  both  with  no  visually  identifiable  interference 
during  robot  motion.  Fiber  optic  sensor  calibration  validates 
the  theoretical  modeling  with  satisfactory  sensing  range  and 
resolution  for  prostate  intervention. 

Keywords:  Optical  Force  Sensor,  Fabry-Perot  Interferometer, 
MRI  Compatibility,  Needle  Driver,  Brachytherapy. 

I.  INTRODUCTION 

Subcutaneous  needle,  catheter  and  electrode  insertion  is 
one  of  the  most  common  minimally  invasive  procedures  [1]. 
Needle  placement  error  can  be  categorized  as  intrinsic  and 
extrinsic  ones.  For  intrinsic  ones,  needle  deflection  due  to 
tissue-needle  interaction  causes  the  deviation  of  needle  tip 
from  the  target.  Intra-  and  post-operative  edema  induces 
implanted  seed  drift  for  procedures  like  brachytherapy.  For 
extrinsic  errors,  perturbations  are  caused  by  patient  move¬ 
ment,  respiratory  motion,  and  external  surgical  tool  caused 
tissue  deformation  (e.g.  ultrasound  probe),  etc.  To  compen¬ 
sate  these  errors  is  one  of  the  major  motivations  of  deploying 
active  needle  steering.  The  proposed  needle  driver  is  capable 
of  steering  bevel  tip  needle  and  active  cannula  while  with  a 
clinical  application  on  prostate  brachytherapy. 

Early  MRI-guided  prostate  robots  focus  on  manual  ac¬ 
tuation.  There  is  active  work  being  developed  in  the  area 

H.  Su,  G.A.  Cole  and  G.S.  Fischer  are  with  Automation  and  Interventional 
Medicine  (AIM)  Robotics  Laboratory,  Department  of  Mechanical  Engineer¬ 
ing,  Worcester  Polytechnic  Institute,  100  Institute  Road,  Worcester,  MA 
01609,  USA  [haosu,  gf  ischer  ]  @wpi  .  edu 

M.  Zervas  and  C.  Furlong  are  with  Center  for  Holographic  Studies  and 
Laser  micro-mechaTronics  (CHSLT)  and  NanoEngineering,  Science  and 
Technology  (NEST),  Department  of  Mechanical  Engineering,  Worcester 
Polytechnic  Institute,  100  Institute  Road,  Worcester,  MA  01609,  USA 


Fig.  1.  (Left)  Physical  prototype  of  6-DOF  piezoelectric  actuated  needle 
placement  robot  consisting  of  needle  driver  module  and  Cartesian  gross 
positioning  module.  The  needle  driver  module  provides  8 cm  insertion 
stroke,  5cm  stylet  retraction  stroke  and  40  revolutions  per  minute  rotation 
speed.  The  Cartesian  gross  positioning  module  provides  8 cm  axial  motion, 
3 cm  elevation  and  4cm  lateral  motion.  (Right)  The  robot  prototype  in  the 
bore  of  a  3T  MRI  scanner  with  a  phantom. 


of  pneumatically  actuated  robotic  devices  [2].  Stoianovici 
et  al.  described  a  MRI-compatible  pneumatic  stepper  motor 
and  applied  it  to  robotic  brachytherapy  seed  placement  [3]. 
Our  previous  work  presented  a  pneumatic  servo  system  and 
sliding  mode  control  [4],  [5].  Kokes  et  al.  [6]  reported  a 
pneumatic  needle  driver  system  for  radio  frequency  ablation 
of  breast  tumors.  Song  et  al.  [7]  reported  a  pneumatically 
actuated  modular  robotic  system  with  parallel  mechanism. 

Pneumatic  actuation  does  have  a  low  level  of  image 
interference,  however  the  scalability,  simplicity,  size  and 
inherent  robustness  of  electromechanical  systems  present  a 
clear  advantage  over  pneumatically  actuated  systems.  To 
this  end,  Chinzei  et  al.  [8]  developed  a  general-purpose 
robotic  assistant  with  ultrasonic  motors.  Goldenberg  et  al. 
[9]  presented  targeting  accuracy  and  MRI  compatibility  tests 
for  a  MRI-guided  robot  employing  ultrasonic  actuators  for 
close-bore  MRI  scanners.  Due  to  unacceptable  signal  noise 
from  the  motor,  the  motor  was  disabled  during  the  scanning. 
Krieger  et  al.  [10]  recently  designed  a  transrectal  prostate 
robot  actuated  by  piezoelectric  motors  with  40%  —  60%  SNR 
reduction  under  motion. 


Brace  finger  against 
hub  to  prevent  motion 


Push  sheath  back 


Seeds 


Seeds  left  in  place 


Fig.  2.  (Top)  Brachytherapy  Needle  from  CP  Medical,  (bottom)  schematic 
of  preloaded  needles:  after  insertion,  the  sheath  is  withdrawn  over  the  stylet, 
leaving  the  seeds  in  the  place  (modified  from  [11]). 
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However,  no  prior  work  has  investigated  piezoelectric 
actuated  robotic  systems  for  needle  steering  under  real¬ 
time  high-field  continuous  MRI.  With  its  many  merits,  in 
a  transperineal  manner  which  alleviates  the  requirement  to 
perform  the  implant  procedures  in  a  different  pose  than  used 
for  preoperative  imaging.  Our  ultimate  overall  goal  is  to 
develop  a  teleoperated  needle  placement  system  consists  of 
a  slave  needle  placement  robot  and  a  fiber  optic  force  sensor 
to  achieve  prostate  intervention  under  continuous  high-field 
MRI.  Hence,  the  contributions  of  the  paper  are  (1)  the 
first  demonstration  of  a  6-DOF  needle  placement  robot  with 
steering  capability  under  real-time  3T  MRI  guidance  with 
less  than  2%  SNR  loss  at  full  speed  during  imaging,  and  (2) 
opto -mechanical  design  of  a  high- re  solution  fiber  optic  force 
sensor  to  measure  needle  insertion  force  and  render  haptic 
display. 

This  paper  is  organized  as  follows:  Section  II  describes 
the  system  requirements  and  mechanism  design  including 
system  architecture,  robot  structure  and  optical  tracking 
frame.  Section  III  presents  the  controller  electrical  design 
and  system  setup.  Fiber  optic  sensing  principle  and  opto¬ 
mechanical  design  are  presented  in  Section  IV.  Phantom 
experiment  in  a  3T  close  MRI  bore  and  sensor  calibration 
are  presented  in  Section  V.  Section  VI  concludes  the  paper 
with  discussion  and  future  work. 

II.  NEEDLE  DRIVER  MECHANISM  DESIGN 

A.  System  Concept  and  Specifications 

Besides  the  MRI  compatibility  constraint,  there  are  fol¬ 
lowing  design  considerations: 

1)  Motion  degree  of  freedom:  3-DOF  motion  needle  driver 
and  3-DOF  Cartesian  gross  positioning  stage  as  shown  in  Fig. 
1.  A  coarse  to  fine  architecture  decouples  the  motion  and 
simplifies  the  kinematics,  while  guaranteeing  high  targeting 
accuracy.  As  shown  in  Fig.  2  top,  the  clinical  18Gauge 
needles  for  prostate  brachytherapy  have  an  inner  stylet  and 
hollow  sheath.  Radioactive  seeds  are  pre-loaded  with  5.5mm 
spacers  between  them  before  staring  the  surgery.  During 
the  insertion,  one  hand  holds  the  cannula  and  the  other 
hand  brace  against  stylet  hub  to  prevent  relative  motion. 
After  insertion,  the  sheath  is  withdrawn  over  the  stylet 
while  leaving  the  seeds  in  place.  To  mimic  the  physician 
preload  needle  type  brachytherapy  procedure,  the  needle 
driver  provides  1-DOF  cannula  rotation  about  its  axis  with  1- 
DOF  translational  insertion.  Another  1-DOF  of  translational 
stylet  motion  is  implemented  to  coordinate  the  motion  with 
respect  to  the  cannula.  The  rotation  motion  of  the  cannula 
may  be  used  for  bevel-based  steering  to  limit  deflection  [12] 
or  may  be  used  for  active  cannula  [13]. 

2)  Operation  in  confined  space:  when  the  patient  lies  in  the 
scanner  bore  with  semilithotomy  position,  the  lateral  space 
between  the  legs  is  around  8cm.  To  fit  into  this  space,  the 
width  of  the  needle  driver  module  has  a  wedge  shape  with 
6cm  front  width  (10cm  long)  and  10cm  back  width  (25cm 
long). 

3)  Sterilization:  only  the  plastic  tip  guide,  collet,  nut  and 
guide  sleeve  have  direct  contact  with  the  needle  and  are 


removable  and  sterilizable. 

B.  System  Architecture 

Three-dimensional  surgical  navigation  software  3D  Sheer 
serves  as  a  user  interface  with  the  robot.  The  navigation 
software  is  running  on  a  Linux-based  workstation  in  the 
scanner’s  console  room.  The  system  workflow  follows  a 
preoperative  planning,  optical  frame  registration,  targeting 
and  verification.  OpenIGTLink  [14]  is  used  to  exchange 
control,  position,  and  image  data.  To  perform  dynamic  global 
registration  between  the  robot  and  scanner,  a  passive  tracking 
the  fiducial  frame  is  integrated  to  the  robot  as  shown  in  Fig. 
3. 

C.  Universal  Needle  Clamping  and  Loading  Mechanism 

To  design  a  needle  driver  that  allows  a  large  variety  of 
standard  needles  to  be  used,  a  new  clamping  device  rigidly 
connect  the  needle  shaft  to  the  driving  motor  mechanism  is 
developed  as  shown  in  Fig.  3.  It  consists  of  three  components 
made  of  ABS  plastic:  collet,  collet  nut  and  collet  screw 
shaft.  This  structure  is  a  collet  mechanism  and  a  hollow 
screw  is  twisted  to  fasten  the  collet  thus  rigidly  locks  the 
needle  shaft  on  the  clamping  device.  The  clamping  device 
is  connected  to  the  rotary  motor  through  a  timing  belt.  An 
eccentric  pulley  tensioner  that  is  concentric  with  the  rotary 
piezoelectric  motor  can  freely  adjust  the  distance  between  the 
motor  and  the  clamping  mechanism.  The  clamping  device 
is  generic  in  the  sense  that  each  collet  can  accommodate  a 
wide  range  of  standard  medical  needle  diameters.  The  overall 
needle  diameter  range  for  three  collets  is  from  25  Gauge 
(0.5144mm)  to  16  Gauge  (1.651mm).  By  this  token,  it  can 
not  only  fasten  brachytherapy  needles  but  also  biopsy  needle 
or  most  other  standard  needles  instead  of  designing  some 
specific  structure  to  hold  the  needle  as  those  in  [15].  The 
plastic  needle  guide  with  quick  release  mechanism,  collet, 
nut  and  guide  sleeve  have  direct  contact  with  the  needle  and 
are  low  lost  and  disposable. 

Fiducial  x7  Optical  encoder 


Fig.  3.  A  exploded  view  of  the  needle  clamping  mechanism,  fiducial 
tracking  frame  and  rotary  motor  fixture  with  timing  belt  tensioner. 

Since  the  linear  motor  for  controlling  the  inner  stylet  with 
respect  to  the  outer  cannula  is  collinear  with  the  collet  and 
shaft  (Fig.  3),  it  is  necessary  to  offset  the  shaft  to  manually 
load  the  needle.  A  brass  spring  preloaded  mechanism  (Fig.  1) 
is  proposed  which  provides  lateral  passive  motion  freedom. 
The  operator  can  pinch  the  mechanism  and  offset  the  top 
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motor  fixture  then  load  needle  and  lock  with  the  needle 
clamping.  This  structure  allows  for  easy,  reliable  and  rapid 
loading  of  standard  needles. 

III.  NEEDLE  DRIVER  ELECTRICAL  DESIGN 

A.  Hardware  Architecture 

The  PiezoMotor  actuators  (Uppsala,  Sweden)  chosen  are 
non-harmonic  piezoelectric  motors,  which  have  two  advan¬ 
tages  over  a  harmonic  drive:  the  noise  caused  by  the  driving 
wave  is  much  easier  to  suppress,  and  the  motion  produced  by 
the  motors  is  generally  at  a  more  desirable  speed  and  torque. 
Optical  encoders  (US  Digital,  Vancouver,  Washington)  have 
been  thoroughly  tested  in  a  3T  MRI  scanner  with  satisfactory 
performance. 

B.  Piezoelectric  Actuator  Driver 

Custom  motor  driver  boards  were  developed  [16],  because 
commercially  available  hardware  to  drive  piezoelectric  mo¬ 
tors  do  not  consider  the  MRI  frequency  interference  problem, 
and  it  is  generally  not  possible  to  drive  the  motors  with 
highly  specific  arbitrary  waveforms  without  interference  to 
the  scanner.  The  driver  is  a  4  channel  high  power  arbitrary 
waveform  generator  designed  to  run  piezoelectric  actuators. 
Waveform  tables  are  loaded  over  USB  or  from  SD  card  by  a 
companion  co-processor  who  is  responsible  for  bootstrapping 
and  provisioning  the  FPGA. 

IV.  FABRY-PEROT  INTERFERENCE  FIBER  OPTIC 
SENSOR 

It  is  reported  in  [17]  that  force  sensing  range  for  prostate 
brachytherapy  is  within  20  Newton  and  a  resolution  of  0.01 
Newton  is  sufficient.  Due  to  the  loss  of  tactile  feedback  in 
a  teleoperated  needle  placement  robot  [18],  a  1-DOF  fiber 
optic  force  sensor  that  measures  in  vivo  needle  insertion 
forces  is  proposed  based  on  our  previous  effort  [19]  to  render 
proprioception  associated  with  brachytherapy. 

A  number  of  fiber  optic  force  sensors  for  MRI  applications 
based  on  light  intensity  modulation  have  been  proposed  [20], 
[21],  to  name  a  few.  Due  to  the  limited  space  in  the  robot 
design,  there  is  a  need  to  miniaturize  the  sensor  while  retain 
the  sensing  range  and  resolution  requirement.  Fiber  Bragg 
Grating  (FBG)  sensors  seem  to  be  a  viable  solution.  FBG 
directly  correlate  the  wavelength  of  light  and  the  change 
in  the  desired  strain.  If  the  fiber  is  strained  from  applied 
loads  then  these  gratings  will  change  accordingly  and  allow 
a  different  wavelength  to  be  reflected  back  from  the  fiber. 
However,  the  costly  optical  source,  FBG  fibers  and  spectral 
analysis  equipment  present  formidable  application  for  medi¬ 
cal  instrumentation.  Fabry-Perot  interference  (FPI)  fiber  optic 
sensor  provides  an  amiable  solution  for  high-resolution  force 
sensing  that  only  relies  on  simple  interference  pattern  based 
voltage  measurement. 

A.  Principle  of  Fabry-Perot  based  Fiber  Optic  Sensor 

In  a  Fabry-Perot  strain  sensor,  light  propagates  through 
a  cavity  containing  semi-reflective  mirrors.  Some  light  and 
transmitted  and  some  is  reflected.  As  shown  in  the  top  of 


Fig.  4.  Fabry-Perot  sensing  principle.  (Top)  light  propagation  in  Fabry- 
Perot  cavity,  (bottom)  resulting  fringe  pattern. 


Fig.  4,  the  distance  between  the  two  fiber  tips  is  generally  on 
the  order  of  nanometers  and,  depending  on  the  gauge  length 
(the  active  sensing  region,  defined  as  the  distance  between 
fusion  welds).  Lcavity  is  the  original  cavity  length.  5  is  the 
change  in  the  cavity  length  from  a  given  load.  The  returning 
light  interferes  resulting  in  black  and  white  bands  known  as 
fringes  (Fig.  4  bottom)  caused  by  destructive  and  constructive 
interference.  The  intensity  of  these  fringes  varies  due  to  a 
change  in  the  optical  path  length  related  to  a  change  in  cavity 
length  when  uni-axial  force  is  applied. 

This  phenomenon  can  be  quantified  through  the  summa¬ 
tion  of  two  waves  [22].  By  multiplying  the  complex  conju¬ 
gate  and  applying  Euler’s  identity,  we  obtain  the  following 
equation  of  reflected  intensity  at  a  given  power  for  planar 
wave  fronts: 

I  —  A\  +  A^  +  2AiA2cos(</>i  —  (p 2 )  (1) 

with  Aiand  A2  representing  the  amplitude  coefficients  of 
the  reflected  signals.  The  above  equation  can  be  changed  to 
represent  only  intensities  by  substituting  A^  =  Ii(i  =  1,2) 
and  —  4>2  =  A0  as 

I  =  h+I2  +  2  V hhcosAcf)  (2) 


LDC:  Laser  Diode  Controller 
PLD  :  Pigtailed  Laser  Diode 
ZT :  Z-axis  Translator 
20XYT  :  20X  Objective  Lens 
BS  :  50:50  Beam  Splitter 
PD :  High-Speed  Photodetector 
DAQ  :  Data  Acquisition  System 
PC  :  Processing  Computer 
FPI :  FPI  Sensor 


Fig.  5.  Schematic  diagram  of  the  opto-mechanical  design  to  implement 
FPI  sensor. 

An  FPI  fiber  optic  strain  sensor  (FISO  Technologies, 
Canada)  was  used  to  evaluate  the  systems  resolution  and 


1585 


19 


potential  integration  into  the  robot.  The  main  component 
of  the  FPI  is  the  sensing  cavity,  measuring  15 .8 /am  wide. 
A  glass  capillary  covering  the  sensing  region  is  fusion 
welded  to  the  fiber  in  two  locations  and  encapsulates  the 
sensor.  There  is  an  air  gap  of  approximately  100.5  /im 
wide.  The  total  length  of  the  FPI  sensor,  including  the  glass 
capillary,  and  bare  fiber  is  approximately  20mm.  Besides 
immune  to  electromagnetic  and  RF  signal  and  substantially 
cheaper  than  FBG,  the  advantages  of  this  sensor  includes: 
1)  static/dynamic  response  capability,  2)  high  sensitivity  and 
resolution,  3)  no  interference  due  to  cable  bending  and  4) 
robust  to  a  large  range  of  temperature  variation  (— 40°~250°) 
due  to  air  gap  insulation  to  the  sensing  region. 


B.  Opto -mechanical  Design 

As  depicted  in  Fig.  5,  the  opto-mechanical  design  of  the 
prototype  begins  with  a  pigtailed  laser  diode  (PLD)  which 
emits  light  in  the  830nm  rand  of  the  infrared  line  with  a 
power  of  lmW.  This  diode  is  controlled  by  a  laser  diode 
controller  (LDC)  (ITC-502,  ThorLabs  Inc,  USA)  which  has  a 
PID  built  in  which  helps  stabilize  the  temperature  and  current 
of  the  diode  when  attached  to  laser  cooler.  The  output  of  the 
pigtailed  laser  that  exits  the  FC  connector  (FC)  at  the  end  of 
the  sensor’s  fiber  is  connected  to  a  Z  axis  translator  (ZT). 
This  Z  axis  translator  helps  focus  the  divergent  light  onto  a 
20X  objective  lens  (Olympus,  Japan)  mounted  to  an  X  —  Y 
axis  translator  (20XYT).  This  collimated  light  is  sent  into  a 
50  :  50  beam  splitter  cube  (BS)  (BS017,  Thorlabs  Inc,  USA) 
where  50%  of  the  light  is  split  towards  the  FPI  sensor  and 
the  other  50%  is  not  used. 

The  light  that  is  sent  to  the  sensor  is  focused  onto  the 
50 p/m  core  of  the  sensor’s  multi-mode  fiber.  This  focusing 
is  accomplished  with  the  help  of  another  20X  objective 
lens  mounted  to  an  X  —  Y  axis  translator  which  focuses 
the  light  onto  the  fiber  core  which  is  able  to  adjust  via  a 
Z-axis  translator  which  has  the  FPI  fiber’s  ST  connector 
(ST)  attached  to  it.  The  light  travels  through  the  fiber  and 
into  the  sensing  cavity  and  then  back  reflects  out  the  same 
optical  axis  it  came  in.  This  back  reflected  light  passes 
through  the  20X  objective  lens  and  is  collimated  into  the 
beam  splitter  and  once  through  the  beam  splitter  the  light 
is  sent  into  the  photodetector  (PD)  (DET10A,  Thorlabs  Inc, 
USA).  The  photodetector’s  output  is  digitized  by  a  16-bit 
data  acquisition  system  (DAQ)  (USB  6229-BNC,  National 
Instruments,  USA)  and  a  processing  computer  (PC)  is  used 
to  calculate  the  strain  values. 

C.  Force  Sensing  Design 

Due  to  negligible  friction  force  between  the  needle  and 
needle  guide,  the  reaction  forces  between  the  mechanism  (top 
plate)  and  the  actuator  drive  rod  is  used  to  measure  needle 
insertion  force  as  shown  in  the  top  of  Fig.  6.  The  beam  to 
hold  the  actuator  rod  has  a  small  1.59mm  groove  to  lay 
the  sensor  that  would  extend  30mm  along  its  side  where 
the  FPI  sensor  could  be  embedded.  The  appropriate  length 
was  provided  to  ensure  that  the  PVC  fiber  covering  would 
be  secured  to  the  top  plate  and  provide  added  durability 


Fig.  6.  (Top)  needle  insertion  force  measurement  based  on  motor 
interaction  force,  (bottom)  finite  element  analysis  of  ABS  top  plate  under 
10  Newton  axial  force. 

to  the  sensor.  Because  each  friction  driven  piezoelectric 
actuator  can  provides  12Newton  force,  the  insertion  trans¬ 
lational  motion  is  provided  by  two  linear  motors.  In  terms 
of  the  insertion  force  range,  lONewton  interaction  force  is 
the  maximum  required  for  each  sensor.  The  finite  element 
analysis  in  Fig.  6  illustrates  the  maximum  strain  lOO^e  under 
10  Newton  axial  force  using  ABS  plastic  material  with  a 
Young’s  Modulus  of  2 GPa  and  Poisson’s  ratio  0.34. 

V.  EXPERIMENTS  AND  RESULTS 

To  demonstrate  the  system  MRI  compatibility  of  this 
architecture  and  the  designed  piezoelectric  driver,  a  series  of 
MRI  phantom  tests  were  performed.  The  sensing  capability 
of  FPI  sensor  was  also  demonstrated  by  experiment. 

A.  MRI  Compatibility  Verification 

The  MRI  compatibility  of  the  needle  placement  robot  was 
demonstrated  in  a  Philips  Achieva  3T  system.  The  phantom 
employed  in  the  experiment  was  a  12cm  diameter  plastic 
tube  filled  with  a  copper  sulfate  solution.  The  motor  and 
encoder  were  placed  immediately  adjacent  to  the  left  side  of 
the  coil.  The  controller  was  placed  approximately  3m  from 
the  scanner  bore. 


Baseline  Motor  On  Difference 


Fig.  7.  Representative  results  showing  the  difference  in  images  obtained 
of  baseline  and  motor  running  conditions.  Different  with  the  results  in  [9] 
and  [10],  this  demonstrates  the  real-time  in  situ  needle  steering  capability. 
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Four  imaging  protocols  were  selected  for  evaluation  of 
compatibility  of  the  system:  1)  diagnostic  imaging  Tl- 
weighted  fast  gradient  echo  (T1  FGE/FFE),  2)  diagnostic 
imaging  T2-weighted  fast  spin  echo  (T2  FSE/TSE),  3)  high¬ 
speed  real-time  imaging  fast  gradient  echo  (FGRE),  and  4) 
functional  imaging  spin  echo-planar  imaging  (SE  EPI).  All 
sequences  were  acquired  with  a  slice  thickness  of  5mm  and 
a  number  of  excitations  (NEX)  of  one.  Three  configurations 
were  evaluated  and  used  in  the  comparison:  1)  baseline 
of  the  phantom  only,  2)  motor  unpowered  with  controllers 
DC  power  supply  turned  on,  3)  motor  on  and  robot  is  in 
motion.  Eight  slices  were  acquired  per  imaging  protocol  for 
each  configuration.  Images  obtained  during  motor  operation 
in  the  scanner  are  subtracted  from  the  baseline  images,  as 
shown  in  Fig.  7.  For  statistical  analysis,  SNR  is  utilized  as 
the  metric  for  evaluating  MRI  compatibility  with  baseline 
phantom  image  comparison  [23].  Statistical  analysis  with 
a  Tukey  Multiple  Comparison  confirms  that  no  pair  shows 
significant  signal  degradation  with  a  95%  confidence  interval. 


t=8s  t=10s 


Fig.  8.  Bevel  tip  needle  insertion  snapshots  during  3T  echo-planar  imaging 
at  0.4  second  interval. 

B.  Needle  Insertion  and  Steering  under  Real-Time  MRI- 
Guidance 

A  series  of  experiments  are  performed  to  evaluate  the  sys¬ 
tem  performance  for  needle  insertion  and  steering  capability 
under  real-time  3T  MRI-guidance. 

The  first  test  is  the  dynamic  needle  insertion.  Gelatin 
(12cm  length,  9 cm  width  and  5cm  thickness)  is  utilized  as  a 
tissue  phantom  for  in  vitro  needle  steering.  The  gelatin  was 
mixed  with  boiling  water  at  a  ratio  of  1  to  1.  A  22 Gauge 
medical  needle  (0.82mm  outer  diameter)  with  45°  bevel  tip 
is  used  for  steering  test.  Functional  imaging  spin  echo-planar 
imaging  (field  of  view  240mm,  echo  time  lms,  repetition 
time  2 ms,  flip  angle  20°)  is  utilized  to  monitor  the  real¬ 
time  needle  motion.  This  imaging  protocol  provides  approx¬ 
imately  2 Hz  update  rates.  Needle  insertion  motion  without 
needle  rotation  is  controlled  by  closed-loop  optical  encoder 
feedback  with  proportional-integral-derivative  controller.  Fig. 
8  depicts  six  bevel  tip  needle  insertion  snapshots  during  3T 


echo-planar  imaging  at  0.4  second  interval.  The  needle  shaft 
and  tip  trajectories  are  clearly  visualized  in  the  phantom 
image  without  major  interference  during  robot  motion. 

In  the  second  test,  the  same  22 Gauge  medical  needle  is 
used  to  demonstrate  the  steering  capability  with  MRI  visual¬ 
ization.  The  bevel  tip  is  rotated  toward  left  before  insertion. 
T2-weighted  fast  spin  echo  (field  of  view  240mm,echo  time 
90ms,  repetition  time  3000ms,  flip  angle  90°)  illustrates  the 
final  needle  shape  and  tip  position.  The  same  procedure  is 
repeated  for  bevel  right  before  insertion  and  the  results  are 
shown  in  Fig.  9.  There  is  no  visually  identifiable  interference 
during  needle  robot  controlled  insertion. 

All  the  three  tests  demonstrate  the  in  situ  piezoelectric 
actuation  capability  in  3T  MRI,  thus  enables  real-time  needle 
steering.  The  compatibility  performance  and  dynamic  needle 
insertion  result  is  significant  comparing  with  the  ones  in 
[9]  and  [10],  which  have  40%  —  60%  SNR  reduction  under 
motion  and  must  interleave  motion  with  imaging. 


Fig.  9.  Bevel  tip  needle  steering  image  in  3T  MRI.  (Left)  bevel  left  needle 
insertion  and  (right)  bevel  right  needle  insertion  . 


C.  Fiber  Optic  Force  Sensor  Calibration 


Calibration  was  performed  by  attaching  the  FPI  to  a 
manufactured  ABS  cantilever  beam.  Strain  on  the  beam  was 
calculated  in  terms  of  the  applied  force  F : 


_  YlFLc 
~  bt3E 


(3) 


where  L  is  the  length  of  the  beam,  c  is  the  distance  from 
the  center  of  the  beam  along  the  ^/-direction,  b  is  the  width 
of  the  base,  t  is  the  thickness,  and  E  is  Young’s  modulus. 

In  order  to  calibrate  the  FPI,  the  relationship  between  the 
intensity  of  light  at  the  output  and  the  strain  was  derived.  A 
hanger  system  was  employed  at  the  end  of  the  cantilever 
beam  to  statically  apply  the  load  in  increments  of  the  5 
grams. 

Recall  in  equation  2,  the  change  in  phase  Acj)  of  the 
intensity  equation  is  equal  to  the  wave  number  ^ ,  multiplied 
by  the  length  of  the  sensing  cavity  region  and  the  strain  in 
the  ^-direction: 


^0  _  xx L cavity)  ^ 

A 

This  value  for  the  change  in  phase  was  substituted  into 
intensity  equation  and  it  is  now  possible  to  predict  the  output 
intensity  of  light  as  a  function  of  the  induced  strain: 


I  =  2/q[1  +  cos( 


(cXx L cavity)  > 


(5) 
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Fig.  10.  Calibration  results  showing  voltage  versus  strain  of  the  FPI  sensor 
together  with  the  theoretical  model. 

The  calibrated  system  can  be  seen  in  the  voltage- strain 
graph  shown  in  Fig.  10.  The  theoretically  predicted  rela¬ 
tionship  is  superimposed  in  the  figure.  The  output  voltage 
follows  a  sinusoidal  pattern  that  repeats  over  an  increasing 
applied  force.  The  discrepancy  between  the  measurement  and 
theoretical  model  is  due  to  the  ambient  light  disturbance  to 
the  opto-mechanical  prototype  which  is  not  shielded  during 
experiment.  A  gage  factor  of  A7A8mv//ie  was  calculated 
and  when  using  a  16  bit  data  acquisition  system. 

VI.  CONCLUSION 

This  paper  presents  the  design  of  a  MRI  compati¬ 
ble  piezoelectric  actuated  6-DOF  robot  integrated  with  a 
high-resolution  fiber  optic  force  sensor  for  image  guided 
brachytherapy.  The  MRI  compatibility  test  of  the  robot  and 
the  calibration  result  of  the  sensor  demonstrate  the  real-time 
piezoelectric  actuation  and  sensing  capability.  The  next  step 
of  this  work  will  focus  on  packaging  the  opto-mechanical 
system  and  attain  robust  and  portable  interface  with  the  robot 
controller.  Needle  steering  with  the  proposed  robot  prototype 
will  be  performed  to  demonstrate  the  targeting  accuracy  in 
live  tissue.  Sensor  hysteresis,  fluctuation  and  drift  would 
be  further  investigated  with  multiple  needle  insertions  and 
retractions. 
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Abstract.  Magnetic  resonance  imaging  (MRI)  can  be  a  very  effective  imaging 
modality  for  live  guidance  during  surgical  procedures.  The  rationale  of  MRI-guided 
surgery  with  robot-assistance  is  to  perform  surgical  interventions  utilizing  ’’real- 
time”  image  feedback  while  minimize  operation  time  and  increase  the  surgical  out¬ 
come.  However,  challenges  arise  from  electromagnetic  compatibility  in  the  high-field 
(1.5T  or  greater)  MRI  environment  and  mechanical  constraints  due  to  the  confined 
close-bore  space.  This  paper  presents  two  MRI-compatible  approaches  for  image- 
guided  transperineal  prostate  needle  placement.  Described  is  the  robotic  mecha¬ 
nism,  actuator  and  sensor  design,  controller  design  and  system  integration  for  a 
pneumatically  actuated  robotic  needle  guide  and  a  piezoelectrically  actuated  needle 
placement  robot.  The  two  degree-of- freedom  (DOF)  pneumatic  robot  with  manual 
needle  insertion  has  a  signal  to  noise  ratio  (SNR)  loss  limited  to  5%  with  alignment 
accuracy  under  servo  pneumatic  control  better  than  0.94mm  per  axis.  While  the 
6-DOF  piezoelectrically  actuated  robot  is  the  first  demonstration  of  a  novel  multi 
piezoelectric  actuator  drive  with  less  than  2%  SNR  loss  for  high-field  MRI  operat¬ 
ing  at  full  speed  during  imaging.  The  preliminary  experiments  in  phantom  studies 
evaluates  system  MRI  compatibility,  workflow,  visualization  and  targeting  accuracy. 


1  Introduction 

Prostate  cancer  is  the  most  common  male  cancer  and  the  second  most  common 
type  of  cancer  in  human.  The  estimated  new  prostate  cancer  cases  (192,280)  in 
2009  account  for  25%  incident  cases  in  men  in  the  United  States  [19].  Each  year 
approximately  1.5  million  core  needle  biopsies  are  performed,  yielding  about  220,000 
new  prostate  cancer  cases.  Over  40,000  brachytherapy  radioactive  seed  implantation 
procedures  are  performed  in  the  United  States  each  year,  and  the  number  is  steadily 
growing.  Prostate  specific  antigen  blood  tests  and  digital  rectal  exams  are  the  two 
major  preliminary  prostate  cancer  diagnosis  methods.  However,  prostate  biopsy  is 
the  conclusive  approach  to  confirm  cancer  diagnosis.  Prostate  brachytherapy  and 
cryotherapy  are  often  used  for  early  therapy. 
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Transrectal  ultrasound  (TRUS)  is  the  current  “gold  standard”  for  guiding  both 
biopsy  and  brachytherapy  due  to  its  real-time  nature,  low  cost  and  ease  of  use.  Fig.  1 
shows  the  traditional  template-based  TRUS-guided  approach  to  brachytherapy  seed 
placement.  However,  TRUS-guided  biopsy  has  a  detection  rate  of  only  20  —  30%  [34]. 
The  drawbacks  are  largely  due  to  the  inherent  limitation  of  ultrasound  imaging  it¬ 
self  and  the  mechanical  template  used  in  the  procedure  to  guide  needles.  Ultrasound 
imaging  is  inferior  to  MRI  for  prostate  cancer  treatment  due  to  its  limited  resolution 
and  inability  to  display  implanted  radiation  seeds.  Furthermore,  ultrasound  probe 
deforms  the  prostate  and  induces  seed  migration.  Since  the  dosimetry  plan  is  usually 
performed  on  the  manual  segmentation  of  pre-operative  ultrasound  images,  the  seed 
placement  accuracy  is  deteriorated  due  to  the  imaging  limitation  and  probe  induced 
tissue  deformation.  On  the  other  hand,  the  template  used  to  guide  needles  in  TRUS 
is  a  6cm  x  6cm  mechanical  grid  with  5mm  x  5mm  space  distance  which  limits  the 
entrance  location  at  the  perineum  and  the  positioning  resolution  of  needle  insertion. 
In  addition,  it  forbids  needle  angulation  to  adjust  needle  orientation  for  better  tar¬ 
geting  or  to  avoid  pubic  arch  interference  (PAI),  thus  restricting  the  eligible  patient 
population. 


Fig.  1.  Traditional  template-based  TRUS-guided  approach  to  prostate  brachyther¬ 
apy  (left  and  middle).  Brachytherapy  needle  from  CP  Medical  and  the  schematic 
procedure  of  preloaded  needles:  after  insertion,  the  sheath  is  withdrawn  over  the 
stylet,  leaving  the  seeds  in  the  place  (modified  from  [37])  (right). 


On  the  other  hand,  the  MRI-based  medical  paradigm  offers  several  advantages 
over  other  imaging  counterparts.  First,  MRI  has  a  couple  of  mechanisms  to  provide 
high-fidelity  soft  tissue  contrast  and  spatial  resolution.  Second,  MRI  is  an  essentially 
three-dimensional  imaging  modality  that  permits  almost  arbitrary  imaging  plane 
selection,  even  in  a  dynamic  manner.  Third,  MRI  produces  non-ionizing  thus  imposes 
no  safety  hazard  to  the  patient  or  practitioner.  The  clinical  efficacy  of  MRI-guided 
prostate  brachytherapy  and  biopsy  was  demonstrated  by  D’Amico  et  al.  at  the 
Brigham  and  Women’s  Hospital  using  a  0.5T  open-MRI  scanner  [8].  MR  images 
were  used  to  plan  and  monitor  transperineal  needle  placement.  The  needles  were 
inserted  manually  using  a  guide  comprising  a  grid  of  holes,  with  the  patient  in  the 
lithotomy  position,  similarly  to  the  TRUS-guided  approach. 
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The  motivation  of  deploying  robotic  system  in  prostate  interventions  comes  from 
three  major  aspects.  First,  robot-assisted  surgery  guarantees  good  geometric  accu¬ 
racy  as  it  relies  upon  rigid  mechanical  structure  and  motion/force  control  of  the 
robot,  thus  overcomes  the  accuracy  limits  of  traditional  mechanical  templates  (5mm) 
and  allows  needle  angulation.  Second,  robots  can  be  designed  with  appropriate  scale 
to  fit  into  open  scanner  bores  or  cylindrical  closed-bores.  Third,  since  robots  are 
stable  and  untiring,  they  can  stay  inside  the  scanner  bore  and  perform  the  interven¬ 
tional  procedures  with  teleoperated  control.  This  eliminates  or  reduces  the  necessity 
of  the  iterative  and  time-consuming  procedure  that  usually  takes  several  cycles  of 
imaging  for  registration  inside  bore,  interventions  out  of  bore  (due  to  space  limit) 
and  confirmation  inside  bore.  In  general,  MRI-guided  surgery  with  robot-assistance 
can  ultimately  minimize  operation  time  while  increases  the  surgical  outcome,  thus 
greatly  reduces  the  equipment  cost  and  overhead. 


(g)  (h)  (i) 


Fig.  2.  MRI  guided  robots  and  assistance  device:  Masamune  et  al.  [22]  (a),  Suther¬ 
land  et  al.  [32]  (b),  Goldenberg  et  al.  [17]  (c),  INNOMOTION  [23]  (d),  Stoianovici 
et  al.  (e),  Chinzei  et.  al  [7]  (f),  Elhawary  et.  al  [11]  (g),  Beyersdorff  et.  al  [3]  (h) 
and  Krieger  et  al.  [21]  (i). 


The  challenges,  however,  arise  from  electromagnetic  compatibility  in  the  high- 
field  (1.5T  or  greater)  MRI  environment  and  mechanical  constraints  due  to  the  con- 
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fined  closed-bore  space.  The  manifestation  of  the  bidirectional  MRI  compatibility 
requires  that  both  the  device  should  not  disturb  the  scanner  function  and  should  not 
create  image  artifacts  and  the  scanner  should  not  disturb  the  device  functionality. 
Thus,  in  the  actuation  level,  conventional  DC/ AC  motors  that  relies  on  electromag¬ 
netism  for  conventional  robotic  actuation  is  not  feasible  in  MRI.  In  the  material 
level,  conductive  materials  such  as  metallic  components  can  induce  heating  thus 
special  means  should  be  taken  to  avoid  the  side-effects  of  these  materials  or  even 
avoid  using  them.  Moreover,  the  confined  physical  space  in  closed-bore  high- field 
MRI  presents  stringent  challenges  for  material  selection  and  mechanical  design. 

Thorough  reviews  of  MRI-compatible  systems  to  date  for  image-guided  interven¬ 
tions  are  presented  by  Tsekos,  et  al.  [36]  and  by  Elhawary,  et  al.  [10].  Fig.  2  shows 
a  collection  of  the  MRI-compatible  assistants  developed  to  date.  Robotic  assistance 
has  been  investigated  for  guiding  instrument  placement  in  MRI,  beginning  with  neu¬ 
rosurgery  [22] .  A  most  recent  MRI-guided  neurosurgical  robot  [32]  was  developed  by 
Sutherland  et.  al.  Goldenberg  et  al.  [17]  proposed  an  MRI-compatible  robotic  system 
for  image-guided  prostatic  interventions.  INNOMOTION  [23]  is  the  first  commer¬ 
cially  available  MRI-guided  robot  that  utilizes  pneumatic  actuation  and  was  recently 
acquired  by  Synthes® .  Stoianovici  et  al.  [28]  described  a  MRI-compatible  pneumatic 
stepper  motor  PneuStep,  which  has  a  very  low  level  of  image  interference.  Chinzei 
et.  al  [7]  developed  a  general-purpose  robotic  assistant  for  open  MRI  that  was  sub¬ 
sequently  adapted  for  transperineal  intraprostatic  needle  placement  [9].  Elhawary 
et.  al  [11]  performed  transrectal  prostate  biopsy.  Beyersdorff  et.  al  [3]  performed  tar¬ 
geted  transrectal  biopsy  in  a  1.5T  MRI  unit  with  a  passive  articulated  needle-guide 
and  have  reported  12  cases  of  biopsy  to  date.  Krieger  et  al.  [21]  presented  a  2-DOF 
passive,  un-encoded,  and  manually  manipulated  mechanical  linkage  to  aim  a  needle 
guide  for  transrectal  prostate  biopsy  with  MRI  guidance.  Song  et.  al  [26]  presented  a 
pneumatic  robot  for  MRI-guided  transperineal  prostate  biopsy  and  brachytherapy. 

Generally,  there  are  four  actuation  principles  for  MRI  applications,  namely  re¬ 
mote  actuation,  hydraulic,  pneumatic  and  ultrasonic/piezoelectric  actuators  [16]. 
A  substantial  armamentarium  of  actuation  methods  has  been  studied  in  a  number 
of  image-guided  interventions  and  rehabilitation  scenarios.  Early  MRI-compatible 
robotic  systems  focus  on  manual  driven  or  ultrasonic  motor  driven  and  the  latter 
cannot  run  during  imaging  due  to  significant  signal  loss.  Remote  actuation  suffers 
from  bulky  structure,  low  bandwidth  and  lower  resolution  and  is  not  preferable  for 
robotic  applications.  Hydraulic  and  pneumatic  actuation  are  considered  as  the  silver 
bullet  for  MRI  applications  in  the  sense  that  these  system  can  completely  get  rid  of 
electrical  and  magnetic  fields  inside  the  scanner  by  using  dielectric  materials,  thus 
providing  high  signal  noise  ratio.  Hydraulic  systems  renders  very  large  output  force, 
but  usually  suffer  from  cavitation  and  fluid  leakage.  Pneumatic  actuation  is  easier  to 
maintain  and  back-drivable  due  to  the  compressibility  of  air  and  is  more  favorable  for 
high  bandwidth  applications,  especially  ideal  for  force  control.  The  disadvantages 
of  pneumatic  actuation  is  also  because  of  compressibility  of  air  and  the  induced 
time  delay,  that  makes  it  difficult  to  control  in  millimeter  level  in  which  the  medi¬ 
cal  robotics  usually  aims  to  control  and  excel  than  human.  Ultrasonic/piezoelectric 
actuation  does  not  rely  upon  magnetism,  but  since  there  is  still  high  frequency  elec¬ 
trical  signal  and  the  piezoelectric  elements  are  often  embedded  inside  ferromagnetic 
or  paramagnetic  materials,  piezoceramic  motors  using  commercially  available  motor 
controllers  have  been  evaluated  by  [11]  which  negatively  impacted  image  quality. 
The  difficulty  arises  from  the  actuator  driving  controller  that  usually  induces  sig- 
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nificant  image  artifact  using  off-the-shelf  control  circuits.  However,  the  scalability, 
simplicity,  size  and  inherent  robustness  of  electromechanical  systems  present  a  clear 
advantage  over  pneumatically  actuated  systems. 

Generally,  no  prior  work  has  considered  the  exploitation  of  modular  robotic  sys¬ 
tems  for  in-situ  needle  placement  in  high-field  closed-bore  MRI,  with  its  many  merits, 
in  a  transperineal  manner  which  alleviates  the  requirement  to  perform  the  implant 
procedures  in  a  different  pose  than  used  for  preoperative  imaging.  This  paper  reviews 
the  design  of  two  approaches  that  comply  with  the  compatibility  while  address  in- 
situ  needle  placement  with  reduced  cost  and  system  complexity.  By  leveraging  the 
pneumatic  and  piezoelectric  actuation  approaches,  the  underlying  philosophy  of  this 
research  effort  aims  to  explore  the  two  distinct  methods  to  demonstrate  two  design 
examples  that  provides  semi-automatic  and  teleoperated  needle  placement  robots 
for  prostate  interventions.  Even  the  focus  of  the  design  targets  prostate  cancer, 
this  interventional  system  is  possible  to  be  adopted  for  other  similar  percutaneous 
applications. 

This  chapter  is  organized  as  follows:  Section  2  analyzes  and  articulates  the  design 
requirement.  Section  3  describes  the  2-DOF  pneumatic  needle  placement  robot  de¬ 
sign  with  manual  insertion  and  Section  4  describes  the  6-DOF  fully  actuated  needle 
placement  robot  design  with  piezoelectric  actuation.  Section  5  presents  the  system 
architecture  and  workflow  utilizing  these  two  robots.  A  number  of  experiments  in 
3  Tesla  closed-bore  scanner  and  analysis  are  elaborated  in  Section  6.  Finally,  a  dis¬ 
cussion  of  the  system  is  presented  in  Section  7  that  summarizes  the  experimental 
insight  and  future  work. 


2  Design  Requirement 

2.1  Surgical  Procedure  Analysis  and  Design  Implication 

Robotic  design  starts  from  the  surgical  procedure  analysis  which  provides  the  guide¬ 
line  and  “optimal”  solution  of  designing  clinically  practical  robots.  This  includes 
the  analysis  of  prostatic  intervention  classification  based  on  interventional  approach 
(transrectal,  transperineal  and  transurinal)  and  the  motion  analysis  of  needle  inser¬ 
tion  during  biopsy,  brachytherapy  and  ablation. 

The  transperineal  approach  is  preferable  for  the  following  reasons  [4].  First, 
transrectal  approach  samples  the  peripheral  zone  (the  most  common  location  of 
the  cancer)  across  its  smallest  dimension  and  has  higher  possibility  of  missing  the 
cancer  than  the  transperineal  approach.  Second,  transrectal  method  does  not  allow 
accessing  some  regions  of  the  prostate  that  are  blocked  by  the  urethra.  Third,  the 
transrectal  technique  is  not  extendible  to  therapy  (brachytherapy  or  cryotherapy) 
that  requires  multiple  needle  placements  because  of  increased  rectal  injury  that 
causes  higher  risk  of  infection. 

As  shown  in  Fig.  1  right,  the  clinical  18Gauge  (1.27mm)  needles  for  prostate 
brachytherapy  have  an  inner  stylet  and  hollow  sheath.  Radioactive  seeds  are  pre- 
loaded  with  5.5mm  spacers  between  them  before  staring  the  surgery.  During  the 
insertion,  one  hand  holds  the  cannula  and  the  other  hand  brace  against  stylet  hub  to 
prevent  relative  motion.  After  insertion,  the  sheath  is  withdrawn  over  the  stylet  while 
leaving  the  seeds  in  place.  To  mimic  the  physician  preload  needle  type  brachyther¬ 
apy  procedure,  the  needle  driver  provides  1-DOF  cannula  rotation  about  its  axis 
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with  1-DOF  translational  insertion.  Another  1-DOF  of  translational  stylet  motion 
is  implemented  to  coordinate  the  motion  with  respect  to  the  cannula.  The  rotation 
motion  of  the  cannula  may  be  used  for  bevel-based  steering  to  limit  deflection  or 
may  be  used  for  active  cannula  [40]. 

2.2  System  Requirement 

The  system’s  principal  function  is  accurate  transperineal  needle  placement  in  the 
prostate  for  diagnosis  and  treatment,  primarily  in  the  form  of  biopsy  and  brachyther- 
apy  seed  placement,  respectively.  The  patient  is  positioned  in  the  supine  position 
with  the  legs  spread  and  raised  as  shown  in  Fig.  3  (left).  The  patient  is  in  a  similar 
configuration  to  that  of  TRUS-guided  brachytherapy,  but  the  typical  MRI  bore’s 
constraint  (<  60cm  diameter)  necessitates  reducing  the  spread  of  the  legs  and  low¬ 
ering  the  knees  into  a  semi-lithotomy  position.  The  robot  operates  in  the  confined 
space  between  the  patient’s  legs  without  interference  with  the  patient,  MRI  scanner 
components,  anesthesia  equipment,  and  auxiliary  equipment  present  as  shown  in  the 
cross-section  shown  in  Fig.  3  (right). 

The  average  size  of  the  prostate  is  50mm  in  the  lateral  direction  by  35mm  in 
the  anterior-posterior  direction  by  40mm  in  length.  The  average  prostate  volume 
is  about  35cc;  by  the  end  of  a  procedure,  this  volume  enlarges  by  as  much  as  25% 
due  to  swelling  [37].  For  our  system,  the  standard  60mm  x  60mm  perineal  window 
of  TRUS-guided  brachytherapy  was  increased  to  100mm  x  100mm,  in  order  to 
accommodate  patient  variability  and  lateral  asymmetries  in  patient  setup.  In  depth, 
the  workspace  extends  to  150mm  superior  of  the  perineal  surface.  Direct  access  to 
all  clinically  relevant  locations  in  the  prostate  is  not  always  possible  with  a  needle 
inserted  purely  along  apex-base  direction  due  to  PAI.  Needle  angulation  in  the 
sagittal  and  coronal  planes  will  enable  procedure  to  be  performed  on  many  of  these 
patients  where  brachytherapy  is  typically  contraindicated  due  to  PAI  as  described 
by  Fu,  et  al.  [14]. 

The  robot  is  situated  upon  a  manual  linear  slide  that  positions  the  robot  in  the 
access  tunnel  and  allows  fast  removal  for  reloading  brachytherapy  needles  or  collect¬ 
ing  harvested  biopsy  tissue.  The  primary  actuated  motions  of  the  robot  include  two 
prismatic  motions  which  replicate  the  DOF  of  a  traditional  template-base  approach, 
and  two  rotational  motions  for  aligning  the  needle  axis  to  help  avoid  PAI  and  critical 
structures.  In  addition  to  these  base  motions,  application-specific  motions  are  also 
required;  these  include  needle  insertion,  canula  retraction  or  biopsy  gun  actuation, 
and  needle  rotation.  The  accuracy  of  the  individual  servo-controlled  joints  is  tar¬ 
geted  to  be  0.1mm,  and  the  needle  placement  accuracy  of  the  robotic  system  itself 
is  targeted  to  be  better  than  1.0mm.  This  target  accuracy  approximates  the  voxel 
size  of  the  MR  images  used  which  represents  the  finest  possible  targeting  precision. 
The  overall  system  accuracy,  however,  is  expected  to  be  somewhat  less  when  effects 
such  as  imaging  resolution,  needle  deflection,  and  tissue  deformation  are  taken  into 
account.  The  MR  image  resolution  is  typically  1mm  and  the  clinically  significant 
target  is  typically  5mm  in  size  (the  precision  of  a  traditional  TRUS  template). 

2.3  MRI  Compatibility  Requirements 

Significant  complexity  is  introduced  when  designing  a  system  operating  inside  the 
bore  of  high- field  1.5-3T  MRI  scanners  since  traditional  mechatronics  materials, 
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B.  System  Requirements 


The  kinematic  requirements  for  the  robot  are  derived  from  the 
workspace  analysis.  A  kinematic  diagram  of  the  proposed  sys¬ 
tem  is  shown  in  Fig.  2.  The  robot  is  situated  upon  a  manual  linear 
slide  that  repeatedly  positions  the  robot  in  the  access  tunnel  and 
allows  fast  removal  for  reloading  brachytherapy  needles  or  col¬ 
lecting  harvested  biopsy  tissue.  The  primary  actuated  motions 
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device  and  a  safety  lock.  The  latter  requires  an  additional  controlled  linear  motion 
to  accommodate  the  cannula  retraction  to  release  the  brachytherapy  seeds.  Rotation 
of  the  needle  about  its  axis  may  be  implemented  to  either  “drill”  the  needle  in  to 
limit  deflect,  or  to  steer  the  needle  using  bevel  steering  techniques  such  as  those 
described  in  [40].  Sterility  has  been  taken  into  consideration  for  the  design  of  the 
end  effectors.  In  particular,  the  portions  of  the  manipulator  and  leg  rest  that  come 
in  direct  contact  with  the  patient  or  needle  will  be  removable  and  made  of  materials 
that  are  suitable  for  sterilization.  The  remainder  of  the  robot  will  be  draped.  An 
alternative  solution  is  to  enclose  the  entire  leg  rest  with  the  robot  in  a  sterile  drape, 
thus  completely  isolating  the  robot  from  the  patient  except  for  the  needle. 


Fig.  4.  Pneumatic  robotic  needle  placement  mechanism  with  two  active  DOF  and 
one  passive,  encoded  needle  insertion.  Dynamic  global  registration  is  achieved  with 
the  attached  tracking  fiducial  (left).  The  configuration  of  the  robot  with  the  con¬ 
troller  in  the  scanner  room  (right). 


3.1  Mechanical  Design 

The  design  of  the  mechanism  is  particulary  important  since  there  are  very  confined 
spaces  and  the  robot  is  constructed  without  the  use  of  metallic  links.  The  design 
was  developed  such  that  the  kinematics  can  be  simplified,  control  can  be  made 
less  complex,  motions  may  be  decoupled,  actuators  can  be  aligned  appropriately, 
and  system  rigidity  can  be  increased.  Based  upon  analysis  of  the  workspace  and 
the  application,  the  following  additional  design  requirements  have  been  adopted: 
1)  prismatic  base  motions  should  be  able  to  be  decoupled  from  angulation  since 
the  majority  of  procedures  will  not  require  the  two  rotational  DOFs;  2)  actuator 
motion  should  be  in  the  axial  direction  (aligned  with  the  scanner’s  axis)  to  maintain 
a  compact  profile;  and  3)  extension  in  both  the  vertical  and  horizontal  planes  should 
be  telescopic  to  minimize  the  working  envelope. 

The  primary  base  DOFs  are  broken  into  two  decoupled  planar  motions.  Motion 
in  the  vertical  plane  includes  100  mm  of  vertical  travel,  and  optionally  up  to  15°  of 
elevation  angle.  This  is  achieved  using  a  modified  version  of  a  scissor  lift  mechanism 
that  is  traditionally  used  for  plane  parallel  motion.  By  coupling  two  such  mecha¬ 
nisms,  as  shown  in  Fig.  5  (left),  2-DOF  motion  can  be  achieved.  Stability  is  increased 
by  using  a  pair  of  such  mechanisms  in  the  rear.  For  purely  prismatic,  both  slides 
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Fig.  3.  This  mechanism  provides  for  motion  in  the  vertical  plane.  Coupling 
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scanner  to  standard  actuators  situated  outside  the  high  field. 
MR-compatible  actuators  such  as  piezoceramic  motors  have 
been  evaluated  in  [15],  [20],  and  [23];  however,  these  are  prone 
to  introducing  noise  into  MR  imaging,  and  therefore,  negatively 


mu  of lVfertical“travel,  ancT optional 

tfeatidepgk  TliHiiiusikWtaiMtili  aem ' 

f(r  i^nt^raiiei 

motion.  By  coupling  two  such  mechanisms,  as  shown  in  Fig.  3, 
2-DOF  motion  can  be  achieved.  Stability  is  increased  by  using  a 
pair  of  such  mechanisms  in  the  rear.  For  purely  prismatic,  both 


teen  demonstrated 


_ ,_T. 

pneumatic)  couplings  [22], 

unnecessary.  Assign  F*.  HSSftdfiSWWeijplfetl 

.  y,  .  .  .  ing  mode  coptrol  (SMC)  with  submillimeter  tracking  accuracy 

only  translational  mpjf^  ^ 

,  W4h?&W6TOf1Mshhfl 

g&9#ri°iif  ^f|r%s#lgMt¥n§Sh^l?^  t^^^^tfeeRftlMdioiaOCfoapiiSJ^arSyseQQPn 

,  .  .  .  ,  .  .  .  .  v.  .  .  because.the  cylinder  bore  is  made  or  glass  and  the  piston  and 

ftfotfea  SffleMlT 

QifsfY;s(fescn6erfs^W(y^vpe0me^arnsi§mtois  &a$4Ms*8diow8s  for  stoilkte^puJiMieMittHy  with 
the  nominarcenfer  position.  Fig.  T shorn  the  mechanism  wnere  „„  n  m  xn  T _ 11o. 

wM'nrite1'  md'ttiUaloAwIsldhefoie  f&frawf^signs,  WfcriHgii 

j$£  SI SMf*ISSEateMe#of ' ste»W|.e«j8B»Of  «*£**«»**  £*f«**w* 

itgsrafeigili.  ffiUtens^tema,  i^KUaderipaijUiadlily 

feiS^1^nPb81iftMB1mtreTil!!^ffW^hl?<diimnateiifaillyrs«afeila)ehigU5yiateetMeaJiy 
insulating,  and  sterilizable  plastic  [e.g., 

u  that  the  tail-sale  state  is  locked  and  applied  airpressure  releases 

C.  Actuator  Design  a  spring-loaded  collet  to  enable  motion.  The  brakes  are  disabled 

The  MRI  environment  places  severe  restrictions  on  the  choice 
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ators.  Many  mechatromc  systems  use  electrodynamic  actuation,  however,  the  very 
nature  of  an  electric  motor  precludes  its  use  in  high-field  magnetic  environments. 
Therefore,  it  is  necessary  to  either  use  actuators  that  are  compatible  with  the  MR  en¬ 
vironment,  or  to  use  a  transmission  to  mechanically  couple  the  manipulator  in  close 
proximity  to  the  scanner  to  standard  actuators  situated  outside  the  high  field.  MR- 
compatible  actuators  such  as  piezoceramic  motors  have  been  evaluated  in  [11,33]; 
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designed  for  the  prostate  intervention  with  the  robot  is  de¬ 
scribed  in  [31].  The  interface  allows  smooth  operation  of  the 
system  throughout  the  clinical  workflow  including  registration, 
planning,  targeting,  monitoring,  and  verification  (see  Fig.  8). 
The  workstation  is  connected  to  the  robot  and  the  scanner's 
console  via  Ethernet.  An  open-source  device  connection  and 
communication  tool  originally  developed  in  virtual  reality  re¬ 
search  [32]  was  used  to  exchange  various  types  of  data  such 
as  control  commands,  positional  data,  and  images  among  the 
components.  Fig.  7  shows  the  configuration  used  in  the  system. 
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top,  and  clamp  down  on  the  rod.  The  design  is  such  that  the  fail-safe  state  is  locked 
and  applied  air  pressure  releases  a  spring-loaded  collet  to  enable  motion.  The  brakes 
are  disabled  when  the  axis  is  being  aligned  and  applied  when  the  needle  is  to  be 
inserted  or  an  emergency  situation  arises. 

Proportional  pressure  regulators  were  the  valves  of  choice  for  this  robot  because 
they  allow  for  direct  control  of  air  pressure,  thus  the  force  applied  by  the  pneumatic 
cylinder.  This  is  an  advantage  because  it  aids  in  controller  design  and  also  has  the 
inherent  safety  of  being  able  to  limit  applied  pressure  to  a  prescribed  amount.  Most 
pneumatic  valves  are  operated  by  a  solenoid  coil;  unfortunately,  as  with  electric 
motors,  the  very  nature  of  a  solenoid  coil  is  a  contraindication  for  its  use  in  an  MR 
environment. With  pneumatic  control,  it  is  essential  to  limit  the  distance  from  the 
valve  to  the  cylinder  on  the  robot;  thus,  it  is  important  to  use  valves  that  are  safe 
and  effective  in  the  MR  environment.  By  placing  the  controller  in  the  scanner  room 
near  the  foot  of  the  bed,  air  tubing  lengths  are  reduced  to  5  m.  The  robot  controller 
uses  piezoelectrically  actuated  proportional  pressure  valves,  thus  permitting  their 
use  near  MRI.  A  pair  of  these  valves  provide  a  differential  pressure  of  ±  100  psi 
on  the  cylinder  piston  for  each  actuated  axis.  A  further  benefit  of  piezoelectrically 
actuated  valves  is  the  rapid  response  time  (4  ms).  Thus,  by  using  piezoelectric  valves, 
the  robot’s  bandwidth  can  be  increased  significantly  by  limiting  tubing  lengths  and 
increasing  controller  update  rate. 


3.3  Robot  Controller  Hardware 

MRI  is  very  sensitive  to  electrical  signals  passing  in  and  out  of  the  scanner  room. 
Electrical  signals  passing  through  the  patch  panel  or  wave  guide  can  act  as  antennas, 
bringing  stray  RF  noise  into  the  scanner  room.  For  that  reason,  and  to  minimize  the 
distance  between  the  valves  and  the  robot,  the  robot  controller  is  placed  inside  of 
the  scanner  room  with  no  external  electrical  connections.  The  controller  comprises 
an  electromagnetic  interference  (EMI)  shielded  enclosure  that  sits  at  the  foot  of  the 
scanner  bed,  as  shown  in  Fig.  4  right;  the  controller  has  proved  to  be  able  to  operate 
3m  from  the  edge  of  the  scanner  bore.  Inside  of  the  enclosure  is  an  embedded  com¬ 
puter  with  analog  I/O  for  interfacing  with  valves  and  pressure  sensors  and  a  field 
programmable  gate  array  (FPGA)  module  for  interfacing  with  joint  encoders  (see 
Fig.  6).  Also,  in  the  enclosure  are  the  piezoelectric  servo  valves,  piezoelectric  brake 
valves,  and  pressure  sensors.  The  distance  between  the  servo  valves  and  the  robot  is 
minimized  to  less  than  5  m,  thus  maximizing  the  bandwidth  of  the  pneumatic  actu¬ 
ators.  Control  software  on  the  embedded  PC  provides  for  low-level  joint  control  and 
an  interface  to  interactive  scripting  and  higher  level  trajectory  planning.  Communi¬ 
cation  between  the  low-level  control  PC  and  the  planning  and  control  workstation 
sitting  in  the  MR  console  room  is  through  a  100-FX  fiber-optic  Ethernet  connec¬ 
tion.  In  the  prototype  system,  power  is  supplied  to  the  controller  from  a  filtered  dc 
power  supply  that  passes  through  the  patch  panel;  a  commercially  available  MR- 
compatible  power  supply  will  be  used  in  future  design  iterations.  No  other  electrical 
connections  pass  out  of  the  scanner  room,  thus  significantly  limiting  the  MR  imaging 
interference. 
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4  Needle  Placement  Robot:  Piezoelectric  Approach 

The  newly  developed  piezoelectric  robot  [31],  shown  in  Fig.  7,  offers  the  capability  of 

real-time  in-situ  needle  steering  in  high-field  MRI.  It  consists  of  a  modular  3-DOF 

needle  driver  with  optical  tracking  frame  coupled  similar  to  that  of  the  pneumatic 

robot.  The  modular  needle  driver  simultaneously  provides  needle  cannula  rotation 

and  independent  cannula  and  stylet  prismatic  motion.  For  experimental  purposes,  it 

is  shown  with  a  generic  MRI-compatible  3-DOF  x  —  y  —  z  stage;  however,  application 

s 
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Fig.  7.  Piezoelectric  robotic  needle  placement  mechanism  with  three  DOF  nee¬ 
dle  rotation  and  collinear  translation  for  prostate  brachytherapy  on  a  three  DOF 
Cartesian  stage(left).  Detailed  design  of  the  needle  driver  platform  (right). 


The  needle  placement  robot  consists  of  a  needle  driver  module  (3-DOF)  and 
Cartesian  positioning  module  (3-DOF).  The  material  is  rapid  prototyped  with  ABS 
and  laser  cut  acrylic.  Considering  the  supine  configuration  and  the  robot  workspace, 
the  width  of  the  robot  is  limited  to  6cm.  As  shown  in  Fig.  7  (left),  the  lower  layer 
of  the  needle  driver  module  is  driven  with  linear  piezoelectric  motor  and  the  upper 
layer  provides  cannula  rotation  motion  and  stylet  prismatic  motion. 

4.1  Mechanical  Design 

To  design  a  needle  driver  that  allows  a  large  variety  of  standard  needles,  a  new 
clamping  device  shown  in  Fig.  7  (right)  rigidly  connects  the  needle  shaft  to  the 
driving  motor  mechanism.  This  structure  is  a  collet  mechanism  and  a  hollow  screw 
made  of  stereolithography  ABS  is  twisted  to  fasten  the  collet  thus  rigidly  locks  the 
needle  shaft  on  the  clamping  device.  The  clamping  device  is  connected  to  the  rotary 
motor  through  a  timing  belt  that  can  be  fastened  by  a  eccentric  belt  tensioner. 
The  clamping  device  is  generic  in  the  sense  that  we  have  designed  3  sets  of  collets 
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and  each  collet  can  accommodate  a  width  range  of  needle  diameters.  The  overall 
needle  diameter  range  is  from  25  Gauge  to  7  Gauge.  By  this  token,  it  can  not 
only  fasten  brachy therapy  needle,  but  also  biopsy  needles  and  most  other  standard 
needles  instead  of  designing  some  specific  structure  to  hold  the  needle  handle. 

Once  a  preloaded  needle  or  biopsy  gun  is  inserted,  the  collet  can  rigidly  clamp 
the  cannula  shaft.  Since  the  linear  motor  is  collinear  with  the  collet  and  shaft,  we 
need  to  offset  the  shaft  to  manually  load  the  needle.  We  designed  a  brass  spring 
preloaded  mechanism  that  provides  lateral  passive  motion  freedom.  The  operator 
can  squeeze  the  mechanism  and  offset  the  top  motor  fixture  then  insert  the  loaded 
needle  through  plain  bearing  housing  and  finally  lock  with  the  needle  clamping. 
This  structure  allows  for  easy,  reliable  and  rapid  loading  and  unloading  of  standard 
needles. 


4.2  Piezoelectric  Actuator  Driver 

The  piezoelectric  actuators  (PiezoMotor,  Uppsala,  Sweden)  chosen  are  non-harmonic 
piezoelectric  motor  which  have  two  advantages  over  a  harmonic  drive:  the  noise 
caused  by  the  driving  wave  is  much  easier  to  suppress,  and  the  motion  produced  by 
the  motors  is  generally  at  a  more  desirable  speed  and  torque.  Even  though  piezoelec¬ 
tric  motor  does  not  generate  magnetic  field,  commercial  motor  driver  boards  usually 
induce  significant  image  artifact  due  to  electrical  noise  according  to  the  most  recent 
result  [20].  A  new  low  noise  driver  was  developed  and  its  architecture  is  shown  in 
Fig.  8 (left)  and  Fig.  8  (right)  shows  the  configuration  of  the  robot  with  the  con¬ 
troller  in  the  scanner  room.  Waveform  tables  are  stored  in  RAM  and  utilized  by  a 


Veoc*rireut 


Fig.  8.  Piezoelectric  actuator  driver  architecture  and  prototype  using  FPGA  gen¬ 
erated  waveform  (left),  and  the  configuration  of  the  robot  with  the  controller  in  the 
scanner  room(right). 
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5  Surgical  Visualization  and  Workflow 

5.1  Surgical  Visualization 

The  user  interface  for  the  robot  is  based  on  3D  Sheer  open-source  surgical  naviga¬ 
tion  software  (3D  Sheer  Software,  http://www.slicer.org).  The  navigation  software 
runs  on  a  Linux-based  workstation  in  the  scanner’s  console  room,  and  communi¬ 
cates  to  the  robot  controller  over  a  fiberoptic  connection.  A  customized  graphical 
user  interface  (GUI)  specially  designed  for  the  prostate  intervention  with  the  robot 
is  described  in  [24].  The  interface  allows  smooth  operation  of  the  system  through¬ 
out  the  clinical  workflow  including  registration,  planning,  targeting,  monitoring  and 
verification  (Fig.  9).  The  workstation  is  connected  to  the  robot  and  the  scanner’s 
console  via  Ethernet.  OpenIGT  Link,  open-source  device  connection  and  communi¬ 
cation  tools  originally  developed  for  virtual  reality  research  [27] ,  are  used  to  exchange 
various  types  of  data  including  control  commands,  position  data,  and  images  among 
the  components.  Fig.  10  shows  the  configuration  used  in  the  system  as  described 
in  [35]. 


Workflow 
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and  loaded  into  the  navigation  software.  Registration  is  performed  between  the 
pre^t){W<Hive  planning  images  and  intra-operative  imaging  using  techniques  such 
as  those  described  by  Haker,  et  al.  [18].  Target  points  for  the  needle  insertion  are 
selected 

mined  target  points  are  selected  in  the  planning  GUI.  Once  the  patient  and  the 

robSff#,?  ,%2^-Ua|ledote 

and  passed  to  the  navigation  software  to  calculate  the  6-DOF  pose  of  the  robot  base 

of  the  phantom  and  targets  selected  as  shown  in  Fig.  3.41.  Target  location  and  Z-frame 
fiducial  location  were  sent  to  the  robot  in  the  scanner’s  RAS  coordinates  (i.e  patient  coor¬ 
dinate  system).  The  robot  then  calculated  the  required  joint  positions.  When  commanded 
to  do  so,  the  robot  aligned  the  needle  with  real-time  feedback  of  the  needle  position  over¬ 
laid  on  the  planning  software  as  shown  in  Fig.  3.42  (left).  The  needle  is  inserted  along  the 
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sent  through  the  network  from  the  navigation  software  to  the  robot  controller.  Af¬ 
ter  the  registration  phase,  the  robot  can  accept  target  coordinates  represented  in 
the  image  (patient)  coordinate  system  in  standard  Right- Anterior- Superior  (RAS) 
coordinates. 

During  the  procedure,  a  target  and  an  entry  point  are  chosen  on  the  navigation 
software,  and  the  robot  is  sent  the  coordinates  and  aligns  the  needle  guide  appro¬ 
priately.  In  the  current  system,  the  needle  is  inserted  manually  while  the  needle 
position  is  monitored  by  an  encoded  needle  guide  and  displayed  in  real-time  on  the 
display.  Needle  advancement  in  the  tissue  is  visualized  on  the  navigation  software  in 
two  comp] 
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Fig.  10.  Diagram  shows  the  connection  and  the  data  flow  among  the  robot  compo¬ 
nents.  OpenIGTLink,  an  open-source  device  communication  tool,  is  used  to  exchange 
control,  position,  and  image  data. 


The  interface  software  enables  “closed-loop”  needle  guidance,  where  the  action 
made  by  the  robot  is  captured  by  the  MR  imaging,  and  immediately  fed  back  to  a 
physician  to  aid  their  decision  for  the  next  action.  The  reason  for  keeping  a  human 
in  the  loop  is  to  increase  the  safety  of  the  needle  insertion,  and  to  allow  for  the 
live  MR  images  to  monitor  progress.  The  robot  fully  aligns  the  needle  as  planned 
before  coming  in  contact  with  the  patient.  If  necessary,  the  placement  is  adjusted 
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responsive  to  the  MR  images.  The  physician  performs  the  insertion  under  real-time 
imaging.  Fig.  9  shows  the  planning  software  with  an  MR  image  of  the  phantom 
loaded  and  real-time  feedback  of  the  robot  position  is  used  to  generate  the  overlaid 
needle  axis  model. 

5.2  Needle  Placement  Robot  Tracking  and  Localization 

Dynamic  global  registration  between  the  robot  and  scanner  is  achieved  by  passive 
tracking  the  fiducial  frame  in  front  of  the  robot  as  shown  in  Fig.  7  (right).  The 
rigid  structure  of  the  fiducial  frame  is  made  of  ABS  and  seven  MR  Spot  fiducials 
(Beekley,  Bristol,  CT)  are  embedded  in  the  frame  to  form  a  Z  shape  passive  fiducial. 
Any  arbitrary  MR  image  slicing  through  all  of  the  rods  provides  the  full  6-DOF  pose 
of  the  frame,  and  thus  the  robot,  with  respect  to  the  scanner.  Thus,  by  locating 
the  fiducial  attached  to  the  robot,  the  transformation  between  patient  coordinates 
(where  planning  is  performed)  and  the  needle  placement  robot  is  known.  To  enhance 
the  system  reliability  and  robust,  multiple  slices  of  fiducial  images  are  used  to  register 
robot  position  using  principal  component  analysis  method.  The  end  effector  location 
is  then  calculated  from  the  kinematics  based  on  the  encoder  positions. 

5.3  System  Workflow 

The  workflow  for  the  system  mimics  that  of  traditional  TRUS-guided  prostate  needle 
insertions,  and  are  as  follows: 

1.  Acquire  pre-procedural  MRI  volume  of  the  patient’s  prostate  and  surrounding 
anatomy. 

2.  Select  specific  needle  tip  targets  as  shown  in  Fig.  9. 

3.  Define  corresponding  needle  trajectories. 

4.  Acquire  an  MR  image  of  the  robot’s  tracking  fiducial. 

5.  Register  the  robotic  system  to  patient/scanner  coordinates. 

6.  Load  the  biopsy  needle  or  pre-loaded  brachytherapy  needle  into  the  robot’s 
needle  driver. 

7.  Send  coordinates  in  patient  coordinates  to  the  robot. 

8.  Automatically  align  the  needle  guide  and  lock  in  place. 

9.  Manually  insert  the  needle  along  prescribed  axis  as  virtual  needle  guide  is 
displayed  on  real-time  MR  images  intersecting  the  needle  axis. 

10.  Confirm  correct  placement. 

11.  Harvest  tissue  or  deliver  therapy. 

12.  Retract  the  needle  guide  and  remove  biopsy  or  brachytherapy  needle  from 
robot. 

13.  Update  surgical  plan  as  necessary  based  on  volumetric  imaging  of  the 
prostate  and  knowledge  of  the  intervention  performed. 

14.  Repeat  for  as  many  needles  as  necessary. 

5.4  Teleoperated  Needle  Insertion 

Manual  insertion  was  the  preferred  technique  due  to  the  need  for  tactile  feedback 
during  the  insertion  phase.  However,  it  was  found  that  the  ergonomics  of  manual  in¬ 
sertion  along  this  guide  proved  very  difficult  in  the  confines  of  the  scanner  bore.  The 
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limited  space  in  closed-bore  high-field  MRI  scanners  requires  a  physical  separation 
between  the  surgeon  and  the  imaged  region  of  the  patient.  To  overcomes  the  loss 
of  needle  tip  proprioception  information,  we  are  developing  a  teleoperated  haptic 
system  with  optical  force-torque  sensor,  to  be  integrated  with  a  3-DOF  robotic  nee¬ 
dle  driver  for  MR- guided  prostate  needle  placement.  The  piezoelectrically  actuated 
robot  incorporates  tactile  feedback  and  teleoperation.  As  shown  in  Fig.  11  ,  an  MRI 
compatible  robot  controller  sits  inside  the  scanner  room  and  communicates  to  the 
navigation  software  and  scanner  interface  running  on  a  laptop  in  the  console  room 
through  a  fiber  optic  connection.  The  optical  force  sensor  interface  is  incorporated 
into  the  in-room  robot  controller  and  the  needle  interaction  forces  are  transmitted 
back  to  the  navigation  software  console  along  with  the  robot  position.  We  integrate 
the  haptic  feedback  device  into  the  navigation  software  framework  to  provide  forces 
to  the  operator  and  control  back  to  the  robot. 

A  direct  force  feedback  architecture  [6]  would  be  used  to  control  the  teleoperated 
needle  placement  system.  In  this  iteration  of  the  system,  we  employ  a  commercially 
available  Novint  Falcon  (Novint  Technologies,  Inc.,  Albuquerque,  NM)  haptic  device 
as  the  master  robot.  It  has  3  position  DOF  and  can  be  used  to  position  the  needle 
in  the  Cartesian  space.  In  the  future  version,  the  haptic  device  would  be  designed 
with  piezoelectric  actuation  and  optical  encoding  as  shown  in  Fig.  11.  It  would  be 
operated  inside  the  scanner  room  for  better  observation  of  patient’s  physiological 
condition  and  other  surgical  procedures.  The  human  operator  positions  obtained 
from  the  haptic  interface  are  used  for  trajectory  generation  and  control  of  the  motion 
of  the  slave  manipulator.  The  slave  robot  in  this  design  is  the  second  generation  of 
the  2-DOF  pneumatically  actuated  robotic  assistant  system  which  overcomes  many 
design  difficulties  and  promises  safe  and  reliable  intraprostatic  needle  placement 
inside  closed  high- field  MRI  scanners.  However,  the  addition  of  force  feedback  allows 
incorporation  of  an  actuated  needle  driver  and  firing  mechanism  and  needle  rotation. 
The  contact  forces  between  needle  and  the  tissue  are  measured  by  the  force/toque 
sensor  [30]  and  further  fed  back  to  the  haptic  device  through  its  interface. 


6  Experiments 

The  first  iteration  of  the  needle  placement  robots  has  been  constructed  and  is  op¬ 
erational.  All  mechanical,  electrical,  communications,  and  software  issues  have  been 
evaluated.  The  current  state  of  the  pneumatic  manipulator  is  two  actuated  DOFs 
(vertical  and  horizontal)  with  an  encoded  passive  needle  insertion  stage.  Evaluation 
of  the  robot  is  in  three  distinct  phases:  1)  evaluation  of  the  MR  compatibility  of 
the  robot;  and  2)  evaluation  of  the  workspace  and  accuracy.  The  piezoelectric  ac¬ 
tuated  robot  is  6-DOF  and  the  MRI  compatibility  has  been  evaluated  using  similar 
protocols.  This  section  reports  the  experimental  setup  and  results  of  the  two  robots. 

6.1  Imaging  Protocols 

Four  imaging  protocols  as  shown  in  Table  1,  were  selected  for  evaluation  of  com¬ 
patibility  of  the  system:  1)  diagnostic  imaging  Tl-weighted  fast  gradient  echo  (T1 
FGE/FFE),  2)  diagnostic  imaging  T2- weighted  fast  spin  echo  (T2  FSE/TSE),  3) 
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Fig.  11.  System  architecture  for  the  master  -  slave  haptic  interface  of  the  piezo¬ 
electric  robot.  The  fiber  optic  force  sensor  and  robot  are  placed  near  the  isocenter 
of  the  MRI  scanner,  the  master  manipulator  is  connected  to  the  navigation  software 
interface,  and  the  two  are  couple  through  the  robot  controller  in  the  scanner  room 
using  a  fiber  optic  network  connection. 


high-speed  real-time  imaging  fast  gradient  echo  (FGRE),  and  4)  functional  imag¬ 
ing  spin  echo-planar  imaging  (SE  EPI).  Details  of  the  scan  protocols  are  shown  in 
Table  1.  All  sequences  were  acquired  with  a  slice  thickness  of  5mm  and  a  number 
of  excitations  (NEX)  of  one.  Three  configurations  were  evaluated  and  used  in  the 
comparison:  1)  baseline  of  the  phantom  only,  2)  motor  powered  with  controllers  DC 
power  supply  turned  on  and  3)  system  servoing  inside  MRI  board.  Three  slices  were 
acquired  per  imaging  protocol  for  each  configuration. 

Compatibility  was  evaluated  on  a  3T  Philips  Achieva  scanner.  A  110  mm,  fluid- 
filled  spherical  MR  phantom  was  placed  in  the  isocenter  and  the  robot  placed  such 
that  the  tip  was  at  a  distance  of  120  mm  from  the  center  of  the  phantom  (a  rep¬ 
resentative  depth  from  perineum  to  prostate),  as  shown  in  Fig.  4(right).  For  the 
pneumatic  robot,  the  phantom  was  imaged  using  four  standard  prostate  imaging 
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Table  1.  Scan  Parameters  for  Compatibility  Evaluation  -  3T  Philips  Achieva 


Protocol 

FOV 

TE 

TR 

Flip  Ang  Bandwidth 

T1W 

240  mm 

2.3  ms 

225  ms 

75° 

751  Hz/pixel 

T2W 

240  mm 

90  ms 

3000  ms 

90° 

158  Hz/pixel 

FGRE 

240  mm 

2.1  ms 

6.4  ms 

50° 

217  Hz/pixel 

EPI 

240  mm 

45  ms 

188  ms 

90° 

745  Hz/pixel 

protocols:  1)  T1W  FFE:  T1  weighted  fast  field  gradient  echo;  2)  T2W  TSE:  T2 
weighted  turbo  spin  echo;  3)  TFE  (FGRE):  “real-time”  turbo  field  gradient  echo. 

A  baseline  scan  with  each  sequence  was  taken  of  the  phantom  with  no  robot 
components  using  round  flex  coils  similar  to  those  often  used  in  prostate  imaging. 
The  following  imaging  series  were  taken  in  each  of  the  following  configurations:  1) 
phantom  only;  2)  controller  in  room  and  powered;  3)  robot  placed  in  scanner  bore;  4) 
robot  electrically  connected  to  controller;  and  5)  robot  moving  during  imaging  (only 
with  T1W  imaging).  For  each  step,  all  three  imaging  sequences  were  performed  and 
both  magnitude  and  phase  images  were  collected. 

The  amount  of  degradation  of  SNR  was  used  as  the  measure  of  negative  effects 
on  image  quality.  SNR  of  the  MR  images  was  defined  as  the  mean  signal  in  a  25  mm 
square  at  the  center  of  the  homogeneous  sphere  divided  by  the  standard  deviation 
of  the  signal  in  that  same  region.  The  SNR  was  normalized  by  the  value  for  the 
baseline  image,  thus  limiting  any  bias  in  choice  of  calculation  technique  or  location. 
SNR  was  evaluated  at  seven  3  mm  thick  slices  (representing  a  25  mm  cube)  at  the 
center  of  the  sphere  for  each  of  the  three  imaging  sequences.  When  the  robot  was 
operational,  the  reduction  in  SNR  of  the  cube  at  the  phantom’s  center  for  these 
pulse  sequences  was  5.5%  for  T1W  FFE,  4.2%  for  T2W  TSE,  and  1.1%  for  TFE 
(FGRE).  Further  qualitative  means  of  evaluating  the  effect  of  the  robot  on  image 
quality  are  obtained  by  examining  prostate  images  taken  both  with  and  without  the 
presence  of  the  robot.  Fig.  12  shows  images  of  the  prostate  of  a  volunteer  placed  in 
the  scanner  bore  on  the  leg  rest.  With  the  robot  operational,  there  is  no  visually 
identifiable  loss  in  image  quality  of  the  prostate. 

6.2  MR  Compatibility  of  Pneumatic  Actuated  Needle  Placement 
Robot 

As  can  be  seen  in  Fig.  12,  the  motors  and  encoders  provide  very  small  visually 
identifiable  interference  with  the  operation  of  the  scanner.  Fig.  14  depicts  one  slice 
of  the  tracking  fiducial  frame  which  provides  the  full  position  information  of  the 
robot.  We  utilize  SNR  as  the  metric  for  evaluating  MR  compatibility  with  baseline 
phantom  image  comparison.  For  comparison,  the  SNR  of  each  configuration  was 
normalized  by  the  average  SNR  of  the  3  baseline  images  for  each  imaging  protocol. 
SNR  of  the  MR  images  is  defined  as  the  mean  signal  in  a  25mm  square  at  the  center 
of  the  homogeneous  sphere  divided  by  the  standard  deviation  of  the  signal  in  that 
same  region  [38] .  The  SNR  was  normalized  by  the  value  for  the  baseline  image;  thus 
limiting  any  bias  in  choice  of  calculation  technique  or  location.  The  technique  used 
for  measuring  SNR  is  equivalent  to  that  described  by  the  NEMA  standard  [2] .  The 
SNR  is  calculated  as  described  in  Method  4  of  the  cited  standard: 
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Fig.  12.  Qualitative  analysis  of  prostate  image  quality  for  pneumatic  robot.  Patient 
is  placed  on  the  leg  support  and  the  pneumatic  robot  sits  inside  of  the  support  tunnel 
inside  the  scanner  bore.  T2  weighted  sagittal  and  transverse  images  of  the  prostate 
taken  when  no  robot  components  were  present  and  when  the  robot  was  active  in 
the  scanner. 


SNR  =  - - -  ,  (1) 

imagenoise 

where,  S  is  the  mean  pixel  value  within  the  Measurement  Region  of  Interest 
(MROI)  and  imagenoise  is  the  standard  deviation  (SD)  of  the  selected  Noise  Mea¬ 
surement  Region  of  Interest  (NMROI) .  The  location  of  the  NMROI  has  a  significant 
impact  on  calculated  SNR;  the  selected  NMROI  in  the  top-left  corner  minimizes 
variance  between  image  slices  and  is  outside  of  ghosting  artifacts  as  described  by 
Firbank,  et  al  [12].  Statistical  analysis  with  a  Tukey  Multiple  Comparison  confirms 
that  no  pair  shows  significant  signal  degradation  with  a  95%  confidence  interval. 

6.3  Accuracy  of  Pneumatic  Actuated  Needle  Placement  Robot 

Accuracy  assessment  is  broken  into  two  parts:  localization  and  placement.  These 
two  must  be  distinguished,  especially  in  many  medical  applications.  In  prostate 
biopsy,  it  is  essential  to  know  exactly  where  a  biopsy  comes  from  in  order  to  be 
able  to  form  a  treatment  plan  if  cancer  is  located.  In  brachy therapy  treatment, 
radioactive  seed  placement  plans  must  be  made  to  avoid  cold  spots  where  irradiation 
is  insufficient;  by  knowing  where  seeds  are  placed,  the  dosimetry  and  treatment  plan 
can  be  interactively  updated.  Based  on  encoder  resolution,  localization  accuracy  of 
the  robot  in  free  space  is  better  than  0.1mm  in  all  directions. 

Positioning  accuracy  is  dependent  on  the  servo  pneumatic  control  system.  The 
current  control  algorithm  for  pneumatic  servo  control  is  based  upon  sliding  mode 
control  techniques.  The  pneumatic  cylinder  was  commanded  to  move  back  and  forth 
between  two  points  and  the  steady  state  error  for  each  move  was  recorded.  The 
average  positioning  accuracy  was  0.26mm  RMS  error.  In  addition  to  point-to-point 
moves,  the  positioning  accuracy  of  the  cylinder  alone  in  free  space  has  been  evaluated 
for  dynamic  trajectory  tracking.  The  cylinder  was  commanded  to  follow  a  0.125Hz, 
20mm  amplitude  sine  wave.  Fig.  13  shows  the  tracking  results  for  the  cylinder. 


42 


zero;  as  shown  in  Fig.  3.17,  there  is  a  lag  in  the  valve  response  when  pressure  increases 

To  appear  in:  Springer  Advances  in  Robotics  and  Virtual  Reality 
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Cylinder  Tracking  Results 
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The  RMS  error  for  this  task  was  0.23mm  which  approximates  the  set  deadband  of 

0.30mi?iJ.2.2  Robotic  System  Accuracy 


P  lace®i<p«fei<Rofo,Qfce  accuracy  results  in  the  previous  section  do  not  directly  translate  into 


The  and 

the  same  metric  (SNR).  As  can  be  seen  in  Fig.  14  (left),  the  motors  and  encoders 
.  i upper  limits  ©f  its  travel  due  to,  both  the- encoder  resolution icqrjrespondina  to  motions  in 
providewery  small  visually  identifiable  interference  with  the  operation  or  the  scanner. 

Fig.  14  (right)  depicts  one  slice  of  the  tracking  fiducial  frame  which  provides  the 
full  position  information  of  the  robot.  Statistical  analysis  with  a  Tukey  Multiple 
Comparison  confirms  that  no  pair  shows  significant  signal  degradation  with  a  95% 
confidence  interval.  This  result  presents  significant  improvement  over  recent  research 
result  in  [20].  It  implies  that  piezoelectric  actuation  is  also  feasible  inside  scanner 
bore  during  imaging,  and  thus  makes  it  a  good  candidate  in  robotic  design. 


7  Discussion 

This  chapter  presents  two  MRI-compatible  approaches  to  image-guided  transper- 
ineal  prostate  needle  placement.  Described  is  the  robotic  mechanism,  actuator  and 
sensor  design,  controller  design  and  system  integration  for  a  pneumatically  actuated 
robotic  needle  guide  and  a  piezoelectrically  actuated  needle  steering  platform.  The 
2-DOF  pneumatic  robot  with  manual  needle  insertion  has  a  SNR  loss  limited  to  5% 
with  alignment  accuracy  under  servo  pneumatic  control  better  than  0.94mm  per 
axis.  While  the  6-DOF  piezoelectrically  actuated  robot  is  the  first  demonstration 
of  a  novel  multi  piezoelectric  actuator  drive  with  less  than  2%  SNR  loss  for  high 
field  MRI  operating  at  full  speed  during  imaging.  The  preliminary  experiments  in 
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Fig.  14.  Representative  results  showing  the  images  obtained  of  baseline  and  system 
servoing  inside  scanner  bore  conditions  (left),  one  slice  of  tracking  fiducial  frame 
besides  a  phantom.  This  result  presents  significant  improvement  over  recent  research 
result  in  [20] (right). 


phantom  studies  evaluates  system  MRI  compatibility,  workflow,  visualization  and 
targeting  accuracy. 

Attaining  an  acceptable  level  of  MR  compatibility  required  significant  exper¬ 
imental  evaluation.  Several  types  of  actuators,  including  piezoelectric  motors  and 
pneumatic  cylinder / valve  pairs,  were  evaluated.  Pneumatic  actuators  have  great  po¬ 
tential  for  MRI-compatible  mechatronic  systems.  Since  no  electronics  are  required, 
they  are  fundamentally  compatible  with  the  MR  environment.  However,  there  are 
several  obstacles  to  overcome.  These  include:  1)  material  compatibility  that  was 
overcome  with  custom  air  cylinders  made  of  glass  with  graphite  pistons;  2)  lack 
of  stiffness  or  instability  that  was  overcome  with  the  development  of  a  pneumatic 
brake  that  locks  the  cylinder’s  rod  during  needle  insertion;  and  3)  difficult  control 
that  was  ameliorated  by  using  high-speed  valves  and  shortening  pneumatic  hose 
lengths  by  designing  an  MRI-compatible  controller.  Pneumatic  actuation  seems  to 
be  an  ideal  solution  for  this  robotic  system,  and  it  allows  the  robot  to  meet  all  of  the 
design  requirements.  Further,  MR  compatibility  of  the  system  including  the  robot 
and  controller  is  excellent  with  no  more  that  a  5%  loss  in  average  SNR  with  the 
robot  operational. 

The  pneumatic  and  piezoelectric  systems  have  been  evaluated  in  a  variety  of 
tests.  The  MR  compatibility  has  shown  to  be  sufficient  for  anatomical  imaging  using 
traditional  prostate  imaging  sequences.  Communications  between  all  of  the  elements, 
including  the  robot,  the  low  level  controller,  the  planning  workstation,  and  the  MR 
scanner  real-time  imaging  interface,  are  in  place.  Initial  phantom  studies  validated 
the  workflow  and  the  ability  to  accurately  localize  the  robot  and  target  a  lesion. 

The  pneumatic  actuation  and  piezoelectric  actuation  presents  complementary 
advantages  and  disadvantages.  Pneumatic  actuation  is  considered  as  the  “ultimate 
solution”  for  MRI  applications  in  the  sense  that  these  system  can  completely  get 
rid  of  electrical  and  magnetic  fields  by  using  dielectric  materials.  However,  due  to 
the  promising  result  implemented  in  the  piezoelectric  robot  and  the  scalability,  sim¬ 
plicity,  size  and  inherent  robustness  of  electromechanical  systems  present  a  clear 
advantage  over  pneumatically  actuated  systems.  The  piezoelectric  actuation  is  con¬ 
sidered  as  the  rule  of  thumb  actuation  method  that  has  been  implemented  as  the 
slave  robot  presented  here.  A  future  development  would  be  augmenting  the  slave 
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robot  to  a  teleoperated  master-slave  system.  The  slave  robot  would  perform  the 
needle  placement  under  the  surgeon  motion  command  while  the  needle  insertion 
force  would  be  measured  using  the  fiber  optic  force  sensors  that  are  being  devel¬ 
oped  [29-31]. 
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Abstract.  Direct  magnetic  resonance  imaging  (MRI)  guidance  during  surgical  in¬ 
tervention  can  provide  many  benefits;  most  significantly,  interventional  MRI  could 
be  used  for  planning,  monitoring  of  tissue  deformation,  real-time  visualization  of 
manipulation,  and  confirmation  of  procedure  success.  Direct  MR  guidance  has  not 
yet  taken  hold  because  it  is  often  confounded  by  a  number  of  issues  including  MRI- 
compatibility  of  existing  surgery  equipment  and  patient  access  in  the  scanner  bore. 
This  paper  presents  a  modular  surgical  system  designed  to  facilitate  the  development 
of  MRI-compatible  intervention  devices.  Deep  brain  stimulation  electrode  placement 
and  prostate  br  achy  therapy  seed  placement  robots  represent  two  examples.  Phantom 
and  human  imaging  experiments  validate  the  capability  of  delineating  anatomical 
structures  in  3T  MRI  during  robot  motion. 


1  Introduction 

Diagnostic  magnetic  resonance  imaging  is  one  of  the  most  effective  imaging  modal¬ 
ities  available  to  the  medical  professional  for  viewing  internal  soft  tissue  structures. 
The  ability  to  use  this  modality  for  live  guidance  during  surgical  procedures  would 
prove  invaluable  for  targeting  or  manipulating  internal  structures  that  are  difficult 
to  reach  for  procedures  including  deep  brain  stimulation  (DBS)  and  percutaneous 
prostatic  intervention. 

Robotic  assistance  for  guiding  instrument  placement  in  MRI  for  neurosurgery 
began  with  Masamune,  et  al.  [9]  and  Chinzei,  et  al.  [2].  Various  methods  are  uti¬ 
lized  to  create  and  control  motion  within  an  MRI  environment,  as  was  outlined  by 
Fischer  et.  al.  [5].  The  authors  have  evaluated  both  pneumatic  [4]  and  piezoelectric 
approaches  [3] .  There  is  additional  work  being  developed  in  the  area  of  pneumatically 
actuated  robotic  devices  such  as  the  PneuStep  [10].  While  this  and  other  pneumatic 
technology  does  have  a  very  low  level  of  image  interference,  the  scalability,  simplicity, 
size  and  inherent  robustness  of  electromechanical  systems  present  a  clear  advantage 
over  pneumatically  actuated  systems.  A  recent  piezoelectric  approach  to  prostatic 
interventions  is  described  by  Krieger,  et  al.  [8],  but  as  with  other  published  piezo¬ 
electric  actuation  schemes,  this  one  causes  unacceptable  image  quality  loss  under 
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motion  and  must  interleave  motion  with  imaging.  While  these  devices  and  others  in 
[1,  6]  do  not  depict  the  entire  state  of  the  art,  they  do  highlight  the  problem  being 
addressed:  the  current  inability  to  generate  and  control  precision  electromechanical 
motion  during  live  high-field  (greater  than  1.5T)  closed-bore  MR  imaging  without  af¬ 
fecting  the  quality  of  the  image.  The  main  contributions  being  discussed  in  this  paper 
are  the  system  architecture  for  the  rapid  development  of  procedure-specific  MRI- 
guided  surgical  intervention  systems,  two  example  implementations  of  this  system 
architecture,  and  experimental  validation  of  their  functionality.  The  two  represen¬ 
tative  procedures  targeted  are  DBS  lead  placement  and  prostate  brachytherapy. 


2  Technical  Approach 

One  of  the  greatest  challenges  to  creating  an  MRI-guided  robotic  system  is  achieving 
the  required  MRI  compatibility  in  terms  of  safety  and  image  quality  [5].  While 
there  are  many  safety  concerns  with  MRI  compatibility  in  terms  of  magnetic  and 
conductive  materials,  employing  material  restrictions  during  the  design  process  can 
easily  and  effectively  minimize  these  risks.  It  has  been  shown  that  a  much  more 
difficult  to  control  compatibility  issue  is  image  interference  [8,  2].  Most  common 
forms  of  electronics,  especially  those  containing  oscillators,  will  generate  noise  in 
the  MR  image  that  will  reduce  the  quality  of  the  images,  sometimes  to  the  point 
that  they  will  not  be  useful  for  therapeutic  medicine. 

Avoiding  the  high  frequency  electrical  noise  problems  are  very  difficult  when 
constructing  an  MRI  compatible  system  with  active  motion.  A  major  difficulty  in 
the  proposition  of  constructing  an  MRI  compatible  actively  driven  system  is  that 
there  are  very  few  “off-the-shelf”  products  that  can  be  used  in  the  sensitive  magnetic 
environment  without  destroying  the  quality  of  images.  The  focus  of  this  research  is 
to  develop  a  series  of  modules  that  will  allow  the  rapid  creation  of  MR  compatible 
actively  driven  systems. 

2.1  System  Architecture 

The  components  of  the  system  presented  are  specifically  constructed  for  use  within 
an  MR  scanner  and  the  system  is  designed  to  be  modular  and  highly  expandable  to 
accommodate  a  wide  variety  of  needs.  As  shown  in  Fig.  1,  the  architecture  is  con¬ 
structed  of  three  primary  components  that  are  independent  of  the  hospital  equip¬ 
ment:  user  workstation,  controller,  and  robotic  mechanism.  The  user  workstation 
(typically  a  laptop  computer)  resides  in  the  MRI  console  room,  while  the  robot  con¬ 
troller  is  in  the  scanner  room  and  the  robotic  device  operates  in  the  MRI  scanner 
bore.  The  system  is  completely  portable  and  not  coupled  to  the  scanner  room  -  the 
only  external  connections  are  a  fiber  optic  cable  that  passes  through  the  waveguide 
and  AC  power  from  the  wall  sockets  in  the  MRI  scanner  room. 

The  user  workstation  resides  in  the  MRI  console  room  and  acts  as  the  primary 
planning  and  navigation  software  interface.  The  system  architecture  supports  a  vari¬ 
ety  of  surgical  planning  platforms  by  supporting  the  OpenIGTLink  communication 
protocol  between  the  workstation  and  the  controller  [14].  Using  the  Sheer  naviga¬ 
tion  software  platform  as  described  in  [13],  MR  images  are  compounded  into  a  three 
dimensional  model  which  is  used  for  procedure  planning  and  interactive  guidance 
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Fig.  1.  Architecture  of  the  modular  MRI-compatible  surgical  system. 

tracking.  The  workstation  communicates  with  the  hospital  imaging  equipment  and 
is  responsible  for  acquiring  the  MR  images  and  optionally  controlling  the  scanner. 

The  navigation  software  sends  target  and  motion  commands  to  the  in-room  robot 
controller,  which  sends  back  pose  and  trajectory  information  about  the  robotic  mech¬ 
anism.  As  configured,  the  navigation  software  need  not  know  the  robot  kinematics 
-  all  commands  are  passed  between  the  planning  workstation  and  robot  controller 
in  patient  (i.e.  RAS)  coordinates.  The  controller  itself  houses  a  majority  of  the  orig¬ 
inal  modules:  motor  drivers,  backplane,  power  supplies,  Faraday  cage,  and  control 
computer.  Typically  the  control  computer  within  the  controller  contains  kinematic 
information  about  the  mechanism  within  the  scanner  bore,  and  interprets  pose  com¬ 
mands  from  the  user  workstation  into  joint  commands  for  the  mechanism.  It  then 
sends  these  position  or  velocity  commands  to  the  motor  drivers,  via  the  backplane. 
The  drivers  then  produce  power  signals  to  operate  actuators  in  the  mechanism,  and 
receives  and  analyzes  encoder  information  for  the  corresponding  degrees  of  free¬ 
dom.  Through  this  system  architecture,  a  varying  plurality  of  degrees  of  freedom 
can  easily  and  quickly  be  integrated  into  a  single  system.  To  increase  modularity, 
the  backplane  has  an  ethernet  interface  that  can  be  directly  coupled  to  the  fiber 
optic  media  converter  (as  opposed  to  the  internal  embedded  control  computer)  to 
alternatively  allow  the  higher  level  kinematic  control  to  be  performed  on  a  computer 
outside  the  scanner  room. 

2.2  System  Modules 

All  of  the  presented  original  system  modules  were  developed  because  the  authors  did 
not  feel  that  any  commercial  off-the-shelf  (COTS)  item  would  accomplish  the  desired 
task  without  causing  signal  degradation  in  the  scanner  images.  One  of  the  more 
common  methods  for  creating  actuated  motion  in  an  MRI  is  the  use  of  piezoelectric 
actuators  due  to  their  inherent  non-magnetic  mode  of  operation,  and  the  high  degree 
of  precision  achievable  by  the  actuators.  Unfortunately,  the  available  drivers  used 
to  operate  these  motors,  while  generally  not  causing  interference  by  their  presence 
in  the  scanner,  do  tend  to  cause  a  large  amount  of  image  degradation  while  the 
motors  they  drive  are  being  driven  [7,  5].  The  following  describes  the  specific  custom 
components  of  the  modular  MRI-compatible  robotic  system,  with  an  example  system 
configuration  shown  in  Fig.  2. 
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Fig.  2.  Example  configuration  of  the  robotic  system  with  the  controller  and  other 
equipment  within  the  scanner  room. 

Piezoelectric  Actuator  Driver  The  primary  objective  during  the  design  of  this  de¬ 
vice  was  to  create  a  generalized  controller  that  could  operate  a  multitude  of  piezo¬ 
electric  actuators.  There  are  many  shapes  and  forms  of  piezoelectric  actuators,  and 
they  generally  fall  into  two  major  categories:  harmonic  and  non-harmonic.  Harmonic 
motors,  also  known  as  ultrasonic  motors,  often  operate  at  higher  speed  and  lower 
torque  than  non-harmonic  motors,  and  utilize  high  voltage  (200-300  volt)  driving 
waveforms  at  high  frequency  (30-50kHz).  Some  examples  of  these  motors  are  the 
Shinsei  ultrasonic  motor  as  described  in  [2],  and  the  Nanomotion  HR  series  of  mo¬ 
tors  described  in  [8].  While  typically  the  drivers  created  for  these  motors  output 
simple  sine  waves,  it  was  desired  to  have  the  ability  to  finely  control  the  shape  of 
the  output  waveform  to  optimize  mechanical  performance  and  MR  compatibility.  In 
order  to  accommodate  the  fine  tuning  of  motor  performance,  the  piezoelectric  driver 
board  was  constructed  with  an  FPGA-based  very  high  speed  (12  mega-samples  per 
second  per  channel)  arbitrary  waveform  generator,  coupled  to  a  high  current  linear 
output  stage,  which  is  operated  with  a  separate  independent  power  rail.  The  high 
degree  of  precision  allows  the  user  to  control  high  frequency  signal  components  be¬ 
ing  passed  through  the  output  stage,  which  could  otherwise  be  expressed  as  noise 
or  artifacts  in  scanner  images. 

While  non-harmonic  motors,  operate  at  a  much  lower  frequency  than  harmonic 
motors  (750Hz  to  3kHz)  they  require  a  high  precision  non-sinusoidal  waveform  to 
be  operated  most  effectively.  Examples  of  these  waveforms  are  shown  in  Fig.  3. 
Because  the  shape  of  the  driving  waveforms  for  non-harmonic  motors  is  not  a  si¬ 
nusoidal  wave,  the  wave  can  naturally  be  broken  down  into  much  higher  frequency 
components.  Due  to  the  possibility  of  these  components  generating  noise  in  scanner 
images,  precision  control  of  these  waveforms  is  coupled  with  low-pass  filtering  on  the 
output.  Combining  this  high  speed  arbitrary  waveform  generator  and  high  current 
output  stage  with  a  microprocessor  capable  of  communicating  via  USB  or  through 
the  Ethernet  backplane  allows  these  driving  units  to  be  rapidly  integrated  into  a 
functional  controller.  This  unit  took  on  the  form  shown  in  Fig.  4. 
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Fig.  3.  Example  waveform  for  non-harmonic  motors  with  desired  input  (left)  and 
Output  from  pi0"7^  Hrivor  Knarrl  /Vi  Wo  rWrOT*  mTl"r»n  +  c  rViia  0000  nf  Waveform 
at  up  to  3kHz 


Velocity  Input 


Fig.  4.  Configuration  of  the  generalized  MR-compatible  piezoelectric  driver  board. 
The  board  performs  all  low-level  control  and  waveform  generation  and  drives  up  to 
four  piezo  elements  in  a  single  or  multiple  motors. 

As  can  be  seen  from  the  diagram,  the  signal  processing  stage  receives  an  input 
in  the  form  of  a  velocity  or  position  setpoint.  This  signal  processing  block  can  then 
create  four  independent  channels  of  driving  waveforms  with  a  data  rate  of  12.5  mega 
samples  per  second  (MSPS)  of  14  bit  resolution.  The  signal  generators  include  a 
high  power  output  amplification  stage,  which  passes  its  signals  out  to  the  actuators 
through  a  7r  filter.  The  boards  also  integrate  position  encoding  (from  differential 
encoder  receivers)  such  that  closed-loop  motion  control  loops  can  take  place  within 
the  FPGA,  distributing  the  real-time  processing  of  the  system. 

Backplane  Signal  Aggregator  As  has  been  referenced  earlier,  a  plurality  of  motor 
drivers  are  used  for  each  of  the  configurations  implemented  with  the  system  architec¬ 
ture.  The  software  structure  of  the  system  uses  a  single  server  to  interface  between 
the  piezoelectric  driver  boards  and  the  control  computer.  A  PIC  microcontroller 
ethernet  device  concatenates  the  lines  of  communication  to  and  from  the  controller 
computer  and  the  driver  boards,  utilizing  a  communication  protocol  known  as  the 
Bowler  Communication  System.  This  structure  is  expressed  through  a  backplane 
connector  containing  a  series  of  board  edge  connecters  to  plug  drivers  into,  an  eth¬ 
ernet  board  connector,  and  output  connectors  for  shielded  cables.  Multiple  config¬ 
urations  of  the  backplane  can  be  customized  for  the  requirements  of  the  particular 
application.  The  current  system  is  based  on  a  10-port  backplane  with  two  68-pin 
shielded,  twisted  pair  VHDCI  cables  connecting  differential  encoder  lines  and  motor 
power  to  a  distribution  block  on  the  robot  (5  axes  each). 
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Power  Supply  A  significant  source  of  noise  when  attempting  to  operate  electrical 
equipment  within  an  MR  scanner  room  can  be  switching  power  converters  supplies. 
One  theory  is  this  noise  is  caused  by  the  high  frequency  operation,  generally  on  the 
order  of  (150kHz  to  600kHz)  of  switching  regulators,  coupled  with  their  rejection 
of  noise  both  on  the  input  and  the  output  of  the  regulators.  Because  of  this,  many 
developers  of  MR  compatible  technologies  utilize  linear  regulators  (which  are  highly 
inefficient),  or  power  supplies  external  to  the  scanner  room  to  create  DC  voltages 
that  are  then  passed  through  the  patch  panel  to  supply  the  in-room  equipment 
(which  requires  scanner-specific  connections  and  an  additional  source  of  noise  in¬ 
filtration).  In  the  proposed  architecture,  a  medical  grade  linear  AC-DC  converter 
generated  48  VDC  from  the  in-room  AC  supply.  To  create  a  power  supply  that 
would  be  capable  of  operating  with  a  high  degree  of  efficiency  (over  85%)  while  still 
maintaining  MRI  compatibility,  a  switching  regulator  was  designed  to  operate  at 
a  switching  frequency  of  approximately  15kHz,  while  still  maintaining  an  output 
ripple  voltage  under  5  mV.  In  addition  to  this,  over  15,000  uF  of  input  capacitance 
was  implemented  in  order  to  minimize  noise  that  could  be  rejected  on  the  input  to 
the  power  supply.  The  design  of  this  power  supply  was  extensively  simulated  under 
varying  loading  conditions  to  ensure  that  it  would  operate  effectively  under  changing 
conditions,  such  as  motors  switching  on  and  off,  or  varying  computational  loads. 

Control  Computer  The  control  computer  in  the  in-room  controller  box  in  general 
is  expressed  as  a  COTS  small  form  factor  computer,  that  has  been  modified  to  be 
powered  off  of  the  supply  rails  generated  by  the  power  supply  stages.  Previous  ver¬ 
sions  of  the  controller  have  used  the  PC  104  form  factor,  the  latest  version  uses  a 
mini-ITX  computer.  In  general,  with  a  grounded  housing  and  continuous  shielded 
case  on  a  computer  such  as  this,  without  the  original  noisy  switching  power  supply, 
these  computers  can  be  operated  within  the  equipment  case  without  much  scanner 
interference  [15].  As  described  earlier,  this  computer  interprets  movement  and  target 
commands  from  the  user  workstation  (using  the  OpenIGTLink)  to  joint  and  actua¬ 
tor  commands  through  the  use  of  kinematic  information  about  the  mechanism  being 
utilized.  In  general,  this  computer  communicates  to  the  backplane  through  an  Eth¬ 
ernet  socket  and  to  the  user  workstation  through  a  fiber  optic  Ethernet  connection. 
To  allow  a  more  distributed  software  architecture,  the  communication  proxy  server, 
trajectory  generator,  and  kinematics  engine  can  be  separate  applications  that  com¬ 
municate  over  a  local  network  socket.  The  fiber  optic  network  connection  between 
the  control  computer  and  user  workstation  is  utilized  to  avoid  passing  any  electrical 
signals  through  the  patch  panel  which  can  lead  to  image  interference. 

EMI  Shielded  Enclosure  As  discussed  before,  signal  integrity  preservation  is  one 
of  the  main  difficulties  facing  the  creation  of  MR  compatible  robotic  systems.  A 
key  step  in  reducing  image  degradation  from  electrical  noise,  is  to  encase  all  elec¬ 
tronics  in  a  continuous  Faraday  cage  to  block  as  much  electromagnetic  interference 
(EMI)  being  emitted  from  the  equipment  as  possible.  This  cage  is  extended  through 
the  shielded  cables  carrying  electrical  signals  out  of  the  cage.  This  enclosure  is  de¬ 
signed  to  support  the  modular  nature  of  the  system  architecture  by  combining  an 
integrated  base  of  support  functionality  such  as  power  rail  generation  and  fiber  op¬ 
tic  communication,  with  a  reconfigurable  equipment  area  as  shown  in  Fig.  5.  By 
making  it  easy  to  change  out  the  equipment  and  system  modules  within  this  case, 
new  devices  can  be  integrated  and  tested  streamlining  the  design  process.  After  the 
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Fig.  5.  Solid  model  of  enclosure  showing  swappable  patch  panels  to  allow  for  differ¬ 
ent  equipment  configurations  (left)  and  a  photograph  of  manufactured  case  (right). 

procedure-specific  system  is  tested  and  verified  as  functional,  a  new  controller,  with 
dedicated  devices  can  be  assembled  with  a  more  compact  construction. 

Tracking  Fiducial  One  of  the  challenges  implementing  image  guided  surgical  pro¬ 
cedures,  is  registering  the  robotic  equipment  with  the  image  space.  To  accomplish 
dynamic  global  registration  between  the  robot  and  the  scanner,  a  passive  tracking 
fiducial  frame  shown  in  Fig.  6  is  implemented  as  described  in  [4].  This  fiducial  frame 
is  constructed  of  ABS  plastic  with  seven  embedded  MR  Spot  fiducials  (Beekley, 
Bristol,  CT).  These  embedded  fiducials  form  a  Z  shape  in  three  different  planes  as 
shown  in  6.  By  utilizing  this  arrangement,  any  arbitrary  MR  image  slicing  through 
all  of  the  rods  can  provide  a  full  6-DOF  pose  of  the  frame,  and  thus  the  robot,  with 
respect  to  the  scanner.  By  locating  the  fiducial  attached  to  the  robot,  the  trans¬ 
formation  between  patient  coordinates  and  the  robot’s  needle  driver  is  known.  By 
transmitting  the  fiducial’s  coordinate  frame  to  the  robot  controller,  the  end  effector 
location  is  then  calculated  from  the  kinematics  based  on  encoder  positions. 


Fig.  6.  Left:  Solid  mode^^^^ckin^fiducial  frame,  Middle:tracking  fiducial  frame 
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>n  is  required  when  placing  the  base 
/tion  in  image  space  with  reference  to 


the  patient.  Because  oktheexpeeted'-repeatedmecessity  for  an  X-Y-Z  type  translation 
to  locate  a  procedure-specific  mechanism  to  a  desired  base  point  with  reference  to  a 
patient,  a  3-axis  translational  stage  was  developed  as  a  system  module.  The  X-Y-Z 
stage  is  a  compact  form  factor  design,  with  the  two  perpendicular  horizontal  degrees 
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of  freedom  being  actuated  by  linear  piezoelectric  motors,  and  the  vertical  stage,  due 
to  the  much  higher  forces  expected  to  be  required,  being  actuated  by  a  scissor-lift 
connected  rotary  motor.  This  is  one  example  of  a  base  for  the  application-specific 
manipulator  that  can  be  customized  for  range  of  motion  and  spatial  constraints  as 
required. 


3  Clinical  Applications:  Two  Example  Procedures 

As  a  show  of  the  systems  modularity,  two  example  interventional  systems  are  be¬ 
ing  developed  in  parallel  by  our  research  group  targeting  two  different  procedures: 
prostate  brachytherapy  and  DBS  electrode  placement.  While  these  procedures  are 
very  different,  because  they  are  both  performed  on  deep,  hard  to  reach  soft  tissue 
structures,  they  could  greatly  benefit  from  MR  guidance  to  improve  performance 
and  reduce  trauma.  Each  application  has  a  custom  end  effector  mechanism  that 
plugs  into  a  common  system  architecture. 


Fig.  7.  Photograph  of  the  remote  center  of  motion  (RCM)  linkage  portion  of  the 
DBS  insertion  mechanism  that  couple  to  the  above-mentioned  translational  motion 
base. 


3.1  DBS  Electrode  Placement  System 

The  most  common  target  for  DBS  electrodes,  the  subthalamic  nucleus(STN),  is  deep 
within  the  brain  and  therefore  cannot  be  directly  seen  during  the  procedure.  In  or¬ 
der  to  accurately  target  the  STN  in  a  traditional  procedure,  a  lengthy  registration 
procedure  is  required  during  which  time  the  patient  is  under  a  considerable  amount 
of  motion.  Coupled  with  the  fact  that  the  MR  images  utilized  to  create  the  surgical 
plan  were  taken  before  the  cranium  is  punctured  and  some  of  the  cerebrospinal  fluid 
is  removed,  there  is  a  clearly  an  opportunity  for  tissue  shift  and  registration  inaccu¬ 
racy.  In  order  to  combat  this  registrational  difficulty,  micro  electrode  recordings  are 
often  used  to  verify  correct  electrode  placement,  though  these  recordings  cannot  be 
made  with  the  DBS  electrodes  themselves.  All  of  these  factors  can  be  reduced  to  the 
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same  problem:  there  is  no  live,  in- situ  confirmation  of  electrode  location  while  the 
electrode  is  being  placed.  Inaccurate  placement  of  electrodes  has  been  shown  to  cause 
side  effects  ranging  from  procedure  ineffectiveness  to  new  psychological  symptoms 
and  even  death  [3] .  By  giving  the  surgeon  access  to  interactively  updated  MR  imag¬ 
ing  during  a  procedure,  there  is  an  opportunity  for  greater  accuracy,  in  addition  to 
greatly  streamlining  and  shortening  the  procedure  by  bypassing  complicated,  multi 
modality  registration  procedures,  and  allowing  the  entire  procedure  to  be  performed 
with  the  patient  remaining  in  one  location.  This  reduces  anesthesia  usage,  trauma 
and  cost  of  the  procedure  as  well.  The  equipment  setup  for  this  proposed  system  is 
shown  in  Fig.  7. 

This  system  was  constructed  utilizing  all  of  the  standard  modules  listed  above, 
with  eight  actuator  drivers  corresponding  to  the  eight  possible  degrees  of  freedom 
(DOF).  In  addition  to  this,  the  mechanism  was  designed  and  constructed  to  be 
kinematically  identical  to  the  Leksell  frame,  the  current  standard  for  DBS  electrode 
insertion,  as  shown  in  Fig.  7.  By  making  the  new  mechanism  kinematically  equivalent 
to  the  currently  used  technology  of  the  Leksell  frame  (a  5-DOF  device),  it  is  the 
authors  hope  to  streamline  the  adoption  of  the  new  system  into  use  by  making  the 
procedural  planning  very  similar  for  the  new  and  old  systems.  Initially,  a  5-DOF 
system  is  targeted  (three  translational  DOF,  two  rotational  DOF),  indicating  a 
manual  insertion,  though  eventually  an  8-DOF  system  will  be  developed  to  integrate 
active  insertion  control,  as  well  as  providing  an  additional  two  degrees  of  rotation 
at  the  end  effector  to  increase  armature  dexterity  as  seen  in  7. 


3.2  Prostate  Needle  Placement  System 

The  needle  placement  mechanism  utilizes  the  same  types  of  actuators  and  encoders 
as  the  DBS  system,  thus  an  identical  controller  configuration  is  utilized.  This  needle 
placement  mechanism,  constructed  primarily  of  rapid  prototyped  ABS  plastic  and 
laser  cut  acrylic,  is  comprised  of  two  main  segments,  a  3-DOF  needle  driving  mod¬ 
ule  (insertion,  cannula  retraction,  and  needle  rotation)  and  the  generalized  3-DOF 
cartesian  positioning  module  as  described  in  [12,  11].  This  system  was  designed  to 
work  in  a  configuration  similar  to  TRUS-guided  brachytherapy,  and  as  such  the  MRI 
bore’s  60  cm  diameter  constrains  the  spread  of  the  legs  and  limited  the  width  of 
the  robot  to  7  cm.  The  lower  motion  platform  layer  provides  linear  motion  from  em¬ 
bedded  piezoelectric  actuators,  while  the  upper  needle  driver  module  layer  provides 
cannula  rotation  and  stylet  prismatic  motion.  This  structure  minimizes  the  “between 
leg”  space  while  the  lower  Cartesian  stage  takes  advantage  of  the  available  “under 
leg”  space. 


4  Experimental  Validation 

Phantom  Test  -  Quantitative  Study  To  date  there  have  been  a  series  of  exper- 
iments  with  this  system,  mainly  focusing  on  the  ability  of  the  system  to  generate 
and  control  motion  while  remaining  MR  compatible.  All  of  the  modules  presented 
earlier  have  been  tested  individually  and  as  a  whole,  to  ensure  that  no  aspect  of 
the  system  creates  an  unacceptable  amount  of  image  interference.  The  specific  scan 
protocols  are  the  same  in  [15].  The  metric  of  choice  to  measure  signal  loss  is  signal 
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Table  1.  Experimental  Results  of  MRI  Compatibility  Evaluation  showing  average 
SNR  for  each  case  and  percent  difference  from  baseline 


Protocol 

Baseline 

Motor  Off  (%change) 

Motor  Running  (%change) 

T1W 

148.7 

150.5  (1.24%) 

149.8  (0.76%) 

T2W 

620.4 

631.8  (1.84%) 

629.4  (1.46%) 

FGRE 

141.2 

142.8  (1.19%) 

141.6  (0.30%) 

EPI 

228.4 

223.6  (2.09%) 

226.3  (0.92%) 

to  noise  ratio  (SNR)  which  compares  an  image  section  that  is  expected  to  be  100% 
signal  with  an  image  section  expected  to  be  100%  dark,  as  can  be  seen  in  Fig.  8. 
The  results  of  the  total  system  demonstrate  a  very  low  amount  of  interference  as 
ca 


Fig.  8.  Left:  Image  subtraction  analysis  of  interference  caused  by  equipment.  Right: 
Plot  of  normalized  SNR  comparing  prototype  presented  in  this  paper  and  the  one 
in  [8]. 


Pre-  Clinical  Evaluation  -  Qualitative  Study  In  another  series  of  tests  performed, 
images  of  a  living  human  subject’s  brain  were  acquired  during  robot  motion  and 
shown  to  several  radiologists;  the  images  were  determined  to  be  indistinguishable 
from  images  taken  without  any  equipment  being  present  in  the  scanner  room.  A 
qualitative  demonstration  of  an  operating  surgical  system  on  image  quality  is  shown 
in  the  anatomy  image  Fig.  9. 
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Fig.  9.  Clear  image  of  patient  anatomy  indicating  usefulness  of  images  for  surgical 
planning.  Note  in  the  image,  the  lack  of  background  noise. 

5  Discussion 

It  has  been  shown  that  the  modular  surgical  system  is  both  capable  of  creating  and 
controlling  motion  in  an  MR  scanner  without  causing  image  interference,  and  that 
the  system  can  be  adapted  to  a  multitude  of  procedure  specific  mechanisms.  Physi¬ 
cal  performance  of  the  system  is  to  be  thoroughly  evaluated  beginning  with  robot’s 
precision  and  clinical  ergonomics  of  the  system’s  configuration.  The  first  phase  of 
this  testing  will  be  to  determine  that  accuracy  with  which  the  system  can  move 
and  control  its  degrees  of  freedom  with  reference  to  itself.  The  next  phase  of  this 
testing  will  be  to  determine  that  accuracy  with  which  the  system  can  position  itself 
with  respect  to  image  space.  Once  these  performance  areas  are  verified,  the  system 
will  be  tested  for  ergonomics  and  usability.  After  the  final  rounds  of  experimental 
validation  discussed  in  the  previous  section  is  completed,  these  two  systems  will 
be  moved  into  cadaver  and  clinical  trials,  while  new  modules  and  procedures  to  be 
incorporated  into  the  modular  device  set  will  be  pursued.  MRI  is  a  highly  effective 
soft  tissue  imaging  system,  and  the  ability  to  utilize  this  procedure  in-vivo  coupled 
with  precision  computer  controlled  motion  will  prove  to  be  an  invaluable  asset  in 
the  future  development  of  minimally  invasive  surgery.  The  authors  hope  the  pre¬ 
sented  modular  robot  controller  system  will  help  to  expedite  the  the  development 
of  clinically  viable  MR  image-guided  robotic  surgery  systems. 

Acknowledgements  This  work  was  funded  in  part  by  the  Congressionally  Directed 
Medical  Research  Programs  Prostate  Cancer  Research  Program  New  Investigator 
Award  W81XWH-  09-1-0191  and  Worcester  Polytechnic  Institute  internal  funds. 

References 

1.  Carpi,  F.,  Khanicheh,  A.,  Mavroidis,  C.,  De  Rossi,  D.:  MRI  compatibil¬ 
ity  of  silicone-made  contractile  dielectric  elastomer  actuators.  Mechatronics, 
IEEE/ASME  Transactions  on  13(3),  370  -374  (2008) 

2.  Chinzei,  K.,  Hata,  N.,  Jolesz,  F.A.,  Kikinis,  R.:  Surgical  assist  robot  for  the 
active  navigation  in  the  intraoperative  MRI:  hardware  design  issues  1,  727-732 
(2000) 


59 


To  appear  in:  Experimental  Robotics  -  Springer  Tracts  in  Advanced  Robotics 


12  G.A.  Cole,  K.  Harrington,  H.  Su,  A.  Camilo,  J.  Pilitsis  and  G.S.  Fischer 

3.  Cole,  G.,  Pilitsis,  J.,  Fischer,  G.S.:  Design  of  a  robotic  system  for  mri-guided 
deep  brain  stimulation  electrode  placement.  In:  Proc.  IEEE  Int.  Conf.  Robotics 
and  Automation  ICRA  ’09,  pp.  4450-4456  (2009) 

4.  Fischer,  G.S.,  Iordachita,  I.,  Csoma,  C.,  Tokuda,  J.,  DiMaio,  S.P.,  Tempany, 
C.M.,  Hata,  N.,  Fichtinger,  G.:  MRI-Compatible  pneumatic  robot  for  transper- 
ineal  prostate  needle  placement.  IEEE/ASME  Transactions  on  Mechatronics 
13(3),  295-305  (2008) 

5.  Fischer,  G.S.,  Krieger,  A.,  Iordachita,  I.,  Csoma,  C.,  Whitcomb,  L.L.,  Fichtinger, 
G.:  MRI  compatibility  of  robot  actuation  techniques-a  comparative  study.  Med 
Image  Comput  Comput  Assist  Interv  11  (Pt  2),  509-517  (2008) 

6.  Gassert,  R.,  Moser,  R.,  Burdet,  E.,  Bleuler,  H.:  MRI/fMRI-compatible  robotic 
system  with  force  feedback  for  interaction  with  human  motion.  Mechatronics, 
IEEE/ASME  Transactions  on  11(2),  216  -224  (2006) 

7.  Kokes,  R.,  Lister,  K.,  Gullapalli,  R.,  Zhang,  B.,  MacMillan,  A.,  Richard,  H., 
Desai,  J.:  Towards  a  teleoperated  needle  driver  robot  with  haptic  feedback  for 
RFA  of  breast  tumors  under  continuous  MRI.  Medical  Image  Analysis  13(3), 
445  -55  (2009) 

8.  Krieger,  A.,  Iordachita,  I.,  Song,  S.E.,  Cho,  N.,  Guion,  P.,  Fichtinger,  G., 
Whitcomb,  L.:  Development  and  preliminary  evaluation  of  an  actuated  MRI- 
compatible  robotic  device  for  MRI-guided  prostate  intervention.  In:  Robotics 
and  Automation  (ICRA),  2010  IEEE  International  Conference  on,  pp.  1066  - 
1073  (2010) 

9.  Masamune,  K.,  Kobayashi,  E.,  Masutani,  Y.,  Suzuki,  M.,  Dohi,  T.,  Iseki,  H., 
Takakura,  K.:  Development  of  an  MRI-compatible  needle  insertion  manipulator 
for  stereotactic  neurosurgery.  J  Image  Guid  Surg  1(4),  242-248  (1995) 

10.  Stoianovici,  D.,  Patriciu,  A.,  Petrisor,  D.,  Mazilu,  D.,  Kavoussi,  L.:  A  new  type 
of  motor:  pneumatic  step  motor.  IEEE/ASME  Transactions  on  Mechatronics 
12(1),  98-106  (2007) 

11.  Su,  H.,  Camilo,  A.,  Cole,  G.,  Hata,  N.,  Tempany,  C.,  Fischer,  G.:  High- field  MRI 
compatible  needle  placement  robot  for  prostate  interventions.  In:  Proceedings  of 
MMVR18  (Medicine  Meets  Virtual  Reality).  Newport  Beach,  California,  USA 
(2011) 

12.  Su,  H.,  Shang,  W.,  Cole,  G.,  Harrington,  K.,  Gregory,  F.S.:  Haptic  system  de¬ 
sign  for  MRI-guided  needle  based  prostate  brachytherapy.  In:  IEEE  Haptics 
Symposium  2010.  IEEE,  Boston,  MA,  USA  (2010) 

13.  Tokuda,  J.,  Fischer,  G.S.,  DiMaio,  S.P.,  Gobbi,  D.G.,  Csoma,  C.,  Mewes,  P.W., 
Fichtinger,  G.,  Tempany,  C.M.,  Hata,  N.:  Integrated  navigation  and  control 
software  system  for  MRI-guided  robotic  prostate  interventions.  Comput  Med 
Imaging  Graph  34(1),  3-8  (2010) 

14.  Tokuda,  J.,  Fischer,  G.S.,  Papademetris,  X.,  Yaniv,  Z.,  Ibanez,  L.,  Cheng,  P., 
Liu,  H.,  Blevins,  J.,  Arata,  J.,  Golby,  A.J.,  Kapur,  T.,  Pieper,  S.,  Burdette, 
E.C.,  Fichtinger,  G.,  Tempany,  C.M.,  Hata,  N.:  Openigtlink:  an  open  network 
protocol  for  image-guided  therapy  environment.  Int  J  Med  Robot  5(4),  423-434 
(2009) 

15.  Wang,  Y.,  Cole,  G.A.,  Su,  H.,  Pilitsis,  J.G.,  Fischer,  G.S.:  MRI  compatibility 
evaluation  of  a  piezoelectric  actuator  system  for  a  neural  interventional  robot. 
In:  Proc.  Annual  Int.  Conf.  of  the  IEEE  Engineering  in  Medicine  and  Biology 
Society  EMBC  2009,  pp.  6072-6075  (2009) 


60 


